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Abstract 


A  new  Aero-Optics  Pressurized  Flow  Facility  was  designed  and  built  to  support  research, 
defence,  and  educational  objectives  in  the  study  of  aero-optical  interactions  in  turbu¬ 
lent  flows.  The  equipment  purchased  as  part  of  this  DURIP  Grant  has  produced  the 
new  aero-optics  Pressurized  Flow  Facility  at  the  University  of  California,  Irvine,  whose 
goals  are:  (a)  to  enable  high-resolution  laboratory  measurements  of  optical  beam  prop¬ 
agation  in  turbulent  flows,  (b)  to  contribute  to  physical  understanding,  modeling,  and 
control  of  aero-optical  flow  phenomena,  (c)  to  contribute  to  the  advancement  of  Air 
Force  aero-optical  applications,  and  (d)  to  educate  engineering  students  in  aeronautical 
applications  of  aero-optics  with  hands-on  laboratory  experience.  The  equipment  selected 
and  purchased  as  part  of  this  laboratory  included:  a  multi-purpose  pressure  vessel  with 
ample  optical  access  to  the  interior  via  five  10-in  diameter  3.25-in  thick  spectrosil  quartz 
windows,  an  air  receiver  that  can  be  pressurized  up  to  3000  psig  needed  to  provide  the 
upstream  driving  pressure  of  the  blow-down  experiments,  flow-management  devices  to 
regulate  the  flow,  a  high-resolution  CCD  camera  used  to  directly  measure  the  flow  field,  a 
pulsed  Nd:YAG  laser  operated  in  the  UV  at  266  nm  with  100  mJ/pulse  at  10  Hz,  a  gaseous 
acetone  seeding  system  used  to  provide  the  acetone  vapor  necessary  for  laser-induced  fluo¬ 
rescence  techniques,  and  a  real-time  data-acquisition  computer.  The  unique  design  of  the 
aero-optics  Pressurized  Flow  Facility  can  produce  novel  high-resolution  measurements  of 
aero-optical  beam  propagation  in  turbulent  flows  at  large  Reynolds  numbers.  These  data 
are  useful  for  the  development  of  accurate  models  that  predict  aero-optical  effects  and 
is  also  useful  to  computational  researchers  by  providing  a  database  documenting  aero- 
optical  effects  through  measurements.  Educationally,  the  requested  equipment  provided 
the  opportunity  for  graduate  and  undergraduate  students  to  acquire  measurement  skills 
useful  for  the  study  of  aero-optical  flow  phenomena. 
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1.  Introduction 


In  order  to  develop  techniques  for  the  prediction  and  control  of  aero-optical  phenomena 
in  turbulent  flows,  an  improved  understanding  is  needed  of  the  large-scale  properties  and, 
for  some  applications,  of  the  small-scale  properties  of  aero-optical  distortions  as  well  as  of 
the  turbulent  refractive  fluid  interfaces  that  generate  these  distortions  (e.g.  [1]  Jumper  & 
Fitzgerald  2001,  and  references  therein).  Practically,  the  large-scale  aero-optical  behav¬ 
ior  is  crucial  in  all  applications  involving  optical  beam  propagation  or  imaging  through 
turbulent  shear  flows,  such  as  the  flows  generated  by  airborne  vehicles.  The  small-scale 
behavior  can  also  be  important  for  those  applications  that  require  high-resolution  and/or 
long-range  optical  imaging  or  beam  propagation.  In  addition,  the  behavior  of  aero-optical 
interactions  at  large  Reynolds  numbers  is  important  and  for  high-speed  flight,  it  is  also 
crucial  to  quantify  Mach-number  effects  on  the  aero-optical  behavior  (e.g.  [2]  Gilbert  & 
Otten  1982,  and  references  therein). 

The  large-scale  organized  flow  behavior  is  known  to  provide  the  dominant  aero-optical 
contributions  in  both  incompressible  turbulent  flows  (e.g.  [3]  Truman  k.  Lee  1990; 

[4]  Wissler  &  Roshko  1992;  [5]  Dimotakis,  Catrakis,  &  Fourguette  2001)  and  weakly- 
compressible  turbulent  flows  (e.g.  [1]  Jumper  &  Fitzgerald  2001).  At  higher  flow  com¬ 
pressibilities,  however,  the  extent  to  which  the  large-scale  flow  behavior  is  organized  is  not 
well  understood  (e.g.  [6]  Smits  &  Dussauge  1996;  [5]  Dimotakis,  Catrakis,  &  Fourguette 
2001).  The  large-scale  properties  of  aero-optical  distortions  at  high  compressibility,  and 
the  manner  in  which  they  are  related  to  the  fluid-interfacial  behavior,  are  also  not  well 
understood. 

Since  the  pioneering  aero-optics  studies  by  Liepmann  (1952  [7],  1952b  [8]),  it  has  been 
recognized  that  one  of  the  central  goals  in  aero-optics  research  is  to  relate  the  optical- 
wavefront  distortions  to  the  fluid-mechanical  behavior.  The  optical-wavefront  distortions 
are  usefully  quantified  by  the  optical  path  length  (OPL),  here  denoted  as  A(x,  t),  for  many 
aero-optics  applications.  The  most  relevant  fluid-mechanical  quantity  in  aero-optics  is 
the  refractive-index  field, 

Co 

c(x,  t) 

where  c(x,  t)  is  the  local  speed  of  light  which  is  always  less  than  or  equal  to  Einstein’s 
universal  speed  of  light  Co  in  vacuum  and  independent  of  the  frame  of  reference.  In  turbu¬ 
lent  flows,  the  refractive-index  field  n(x,  t)  can  be  highly  nonuniform  particularly  at  large 
Reynolds  numbers  and  at  high  compressibility.  The  propagation  of  optical  wavefronts 
through  refractive-index  fields  is  governed  by  the  eikonal  equation, 

I V  A|  =  n,  (2) 

where,  as  stated  above,  A  denotes  the  OPL.  The  eikonal  equation  is  useful  to  describe 
fluid-optical  interactions  as  long  as  the  wavelength  of  light  is  smaller  than  the  smallest 
fluid-mechanical  (turbulent)  scale  and  as  long  as  the  optical-beam  energies  are  low  enough 
to  not  change  the  local  refractive  index.  An  experimental  flow  facility  that  can  enable 
the  simultaneous  measurements  of  turbulent  flows  and  flow-generated  optical  wavefront 
distortions  is  crucial  for  understanding  the  aero-optical  interactions.  Relating  the  optical- 
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wavefront  distortions  to  the  fluid-mechanical  behavior  is  one  of  the  major  goals  of  the 
aero-optics  Pressurized  Flow  Facility  built  in  part  by  this  DURIP  Grant. 

2.  Description  of  the  New  Aero-Optics  Pressurized  Flow  Facility 

The  new  Aero-Optics  Pressurized  Flow  Facility  was  designed  and  built  to  allow  direct 
simultaneous  measurements  of  the  flow  field  and  the  flow-generated  optical  wavefront  dis¬ 
tortions.  The  ability  to  simultaneously  measure  the  flow  field  and  the  optical  wavefronts 
distortions  is  one  of  the  most  important  contributions  of  the  new  Aero-Optics  Pressurized 
Flow  Facility  which  is  housed  in  the  Aeronautics  and  Fluid  Dynamics  Laboratories  at 
the  UC  Irvine.  The  new  facility  shown  in  the  photograph  of  figure  1  is  a  pressurized  wind 
tunnel  operated  in  a  blow-down  configuration.  Figure  2  shows  a  drawing  of  the  driver 
reservoir  vessel,  and  figure  3  shows  a  drawing  of  the  main  Aero-Optics  vessel. 


Figure  1:  The  aero-optics  Pressure  Vessel  at  UC  Irvine  designed  for  simultaneous  imag¬ 
ing  of  the  flow  field  and  the  optical  wavefront  distortions.  The  facility  can  generate 
aero-optical  flows  at  large  Reynolds  numbers  Re  ~  106  —  108  and  over  a  wide  range  of 
compressibilities  0.15  <  Mc  <  1.0. 

An  air  reservoir  with  an  internal  volume  of  200-gallon  was  designed  according  to  the 
ASME  Pressure  Vessel  Code  Section  VIII  Division  I  and  it  was  issued  a  U-stamp  upon 
completion  of  it’s  construction  and  inspection.  A  drawing  of  the  air  reservoir  is  shown 
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in  figure  3.  The  air  reservoir  has  a  maximum  pressure  rating  of  3000  psig  in  order  to 
provide  the  necessary  upstream  driving  pressure  to  run  experiments  at  high  Reynolds 
numbers  (Re  ~  106  -  108)  and  across  a  wide  range  of  compressibilities  (Mc  ~  0.15  — 1.0). 
Airflow  from  the  air  reservoir  to  the  experiments  is  delivered  through  high-pressure  6-inch 
diameter  piping.  The  facility  is  operated  in  a  blowdown  configuration  and  the  mass  flow 
rate  to  the  experiments  is  controlled  using  a  pneumatically-actuated  fast-acting  valve  in 
combination  with  fixed-area  metering  devices. 


SIDE  VIEW 


Figure  2:  Drawing  of  the  air  reservoir  which  has  an  internal  volume  of  200  gallons  and 
can  be  pressurized  up  to  3000  psig. 

The  airflow  from  the  upstream  high-pressure  air  reservoir  is  received  by  the  aero-optics 
pressure  vessel  which  has  a  4-foot  inside  diameter  and  8-foot  internal  height.  There  are 
two  6-in  diameter  inlets  located  on  the  bottom  of  the  aero-optics  vessel  and  four  2-in 
outlets  located  on  the  top.  This  aero-optics  vessel  provides  a  controlled  environment  for 
the  experiment’s  test  section  with  the  capability  to  vary  the  internal  pressure  up  to  a 
maximum  of  325  psig  or  20  atm. 
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Figure  3:  Drawing  of  the  main  Aero-Optics  vessel  with  8-foot  internal  height  and  4-foot 
inside  diameter,  with  the  quick-opening  2-foot  diameter  door. 

The  variable  pressure  capability  allows  for  elevated  pressures  in  the  experiment  test- 
section,  which  enhances  the  measurable  aero-optical  effect  and  the  signal-to-noise  ratio 
of  the  flow  and  wavefront  imaging  techniques.  A  2-foot  diameter  custom-designed  breach 
lock  type  door  that  is  hydraulically  operated  permits  access  to  the  interior  of  the  vessel 
for  the  purpose  of  installing  and  adjusting  the  components  of  the  experiment  flow  tunnel. 
The  door  on  the  side  of  the  aero-optics  vessel  is  equipped  with  an  electronic  start/stop 
station  that  enables  the  opening/closing  of  the  door  within  10  seconds.  In  addition,  the 
door  includes  a  safety  inte'rlock  to  prevent  the  opening  of  the  door  when  the  vessel  is 
pressurized.  Optical  access  to  the  interior  of  the  aero-optics  vessel  is  enabled  through 
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five  10.00+/-0.010  inch  diameter  3.25+/-0.010  inch  thick  windows.  The  material  of  the 
windows  was  chosen  to  be  spectrosil  quartz  in  order  to  allow  more  than  90%  transmission 
of  ultraviolet  laser  beams.  The  spectrosil  quartz  windows  were  mounted  on  the  pressure 
vessel  using  custom-designed  weld  units  that  use  compression  screws  to  maintain  the 
windows  in  place  under  the  maximum  pressure  rating  of  325  psig.  The  spectrosil  quartz 
windows  are  polished  on  both  sides  to  a  40/20  surface  finish  and  a  flatness  of  less  than 
1/4  waves  per  inch  over  85%  of  clear  aperture. 

The  aero-optics  Pressure  Vessel  provides  a  controlled  variable  pressure  environment  for 
the  test-section  of  the  experiment.  The  relatively  large  size  of  the  aero-optics  Pressure 
Vessel  enables  the  flow  tunnel  internals  to  be  modified  in  order  to  accommodate  different 
flow  geometries,  e.g.  dual-stream  shear-layer,  single-stream  shear-layer  in  planar  and 
curved  configurations.  As  part  of  this  grant  a  flow  tunnel  of  the  single-stream  shear  layer 
was  designed  and  built.  The  configuration  of  the  single-stream  shear  layer  was  designed 
to  generate  a  subsonic  single-stream  shear-layer  with  a  splitter-plate  test-section  exit 
height  of  3  inches  and  a  span  of  12  inches  as  shown  in  figure  4.  This  configuration 
allows  for  a  comprehensive  study  of  a  separated  flow  of  practical  relevance  for  aero-optics 
where  the  effects  of  Mach  number,  Reynolds  number,  and  density-ratio  may  be  isolated. 
High-Reynolds  number  (Re  ~  106  —  108)  and  high-Mach  number  ( Mc  ~  1.0)  flows  are 
attainable  with  the  current  single-stream  shear-layer  configuration. 


Figure  4:  Schematic  of  the  single-stream  shear-layer  tunnel  insert  used  to  generate  the 
separated  turbulent  shear  layer,  shown  in  the  subsonic  configuration. 

A  single-stream  shear-layer  tunnel  with  a  test-section  entrance  area  of  3-inch  x  12-inch  is 
contained  inside  the  main  pressurized  vessel.  The  configurations  of  the  single-stream  and 
dual-stream  shear  layers  share  many  design  aspects.  The  current  configuration  of  the 
single-stream  shear  layer  consists  of  a  12-inch  x  12-inch  flow  management  section,  which 
contains  two  sets  of  metal  screens  and  two  sets  of  honeycomb  sections.  For  subsonic 
experiments,  after  the  honeycomb  sections  a  converging  nozzle  leads  the  flow  to  the 
test  section  entrance  of  3-inch  height  x  12-inch  span.  To  meet  the  flow  requirements 
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of  supersonic  experiments,  the  flow  management  section  and  contraction  section  are 
adaptable  to  geometrically  smaller  cross-sections  that  also  include  a  converging-diverging 
nozzle.  The  test  section  consists  of  two  fixed  sidewalls  and  an  adjustable  guidewall  that 
spans  the  2-foot  length  of  the  test-section.  The  fourth  plane  of  the  test  section  is  open 
without  restrictions  in  order  to  allow  entrainment  of  the  surrounding  ambient  air.  As 
shown  in  figure  3,  an  entrainment  wall  is  positioned  near  the  test-section  entrance  in 
order  to  provide  near-field  boundary  conditions  consistent  with  previous  studies.  The 
test-section  boundary  conditions  of  the  single-stream  shear  layer  provide  a  favorable 
geometry  for  the  maximum  growth  of  separated  layers  that  begin  to  form  at  the  entrance 
of  the  test-section.  The  maximum  growth  rate  is  a  favorable  characteristic  of  separated 
shear  layers  in  terms  of  capturing  the  flow  structure  using  digital  imaging  techniques. 

A  pulsed  Nd:YAG  laser  provides  the  illumination  source  needed  for  the  measurement  of 
the  flow  images  and  the  optical  wavefront  distortions.  The  Nd:YAG  laser  is  a  Spectra- 
Physics  INDI  model  equipped  with  a  fourth-harmonic  generator  that  enables  its  operation 
in  the  ultraviolet  range  at  266- nm.  The  energy  output  is  /  100  mJ/pulse  at  a  repetition 
rate  of  10  Hertz.  High-quality  optical  mirrors  and  laser-sheet-forming  optics  with  anti- 
reflective  coatings  are  placed  externally  to  the  laser  to  guide  the  optical  beam  to  the 
experiment’s  test  section  trough  the  spectrosil  quartz  windows. 

3.  Facility  Operation  for  Quantitative  Aero-Optical  Imaging 

The  flow  facility  and  shear  layer  tunnel  configuration  described  above  are  used  to  generate 
a  separated  flow  with  the  maximum  growth  rate.  The  optical  wavefront  distortions  can  be 
measured  directly  by  employing  a  Shack-Hartmann  wavefront  sensor.  Direct  quantitative 
measurements  of  the  density  field  may  be  achieved  using  laser-imaging  techniques  such 
as  planar  laser-induced  fluorescence  (PLIF)  of  acetone  vapor.  However,  the  relatively 
high  signal  obtained  from  laser  fluorescence  of  acetone  vapor  compared  to  other  imaging 
techniques  makes  it  a  prospective  candidate  for  flow  imaging  of  separated  flows  (Lozano, 
Yip,  &  Hanson  1992  [9]).  A  density-field  imaging  technique  is  proposed  using  laser- 
induced  fluorescence  of  acetone  where  calibrated  concentrations  of  seeded  acetone  vapor 
are  introduced  in  the  freestream  and  in  the  test-section  ambient  air.  The  calibrated 
concentrations  of  the  two  streams  allows  for  direct  measurements  of  intensity  variations 
due  to  compressibility  effects  without  the  influence  of  molecular  mixing.  This  density 
imaging  technique  provides  a  means  for  separating  the  effects  of  molecular  mixing  of 
dissimilar  gases  from  the  effects  observed  due  to  air  compressibility  at  high  speeds. 

The  objective  of  this  flow  imaging  technique  is  to  directly  relate  the  captured  flow  images 
of  acetone  fluorescence  to  the  index-of-refraction  flow  field.  The  method  for  obtaining 
the  density  flow  field  consists  of  forming  the  output  beam  of  a  Nd:YAG  pulsed  laser  into 
a  thin  two-dimensional  laser  sheet.  The  Nd:YAG  laser  is  operated  in  UV  mode  at  a 
wavelength  of  266-nm  in  order  to  remain  within  the  absorption  spectrum  needed  for  the 
fluorescence  of  acetone  vapor.  The  UV  laser  was  operated  in  10-Hz  mode  with  typical 
energy  output  of  100-mJ  per  pulse.  The  concentration  of  acetone  vapor  was  introduced 
by  means  of  a  metered  bubbling  process  separately  to  the  upstream  air  reservoir  and 
the  test  section  that  is  housed  inside  the  aero-optics  vessel.  This  acetone  vapor  seeding 
procedure  allows  for  the  acetone  vapor  concentration  to  be  matched  in  the  flowing  stream 
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and  in  the  experiment  test  section.  This  acetone  seeding  technique  allows  acetone  vapor 
seeding  of  only  one  of  the  streams,  which  at  weakly  to  highly-compressible  flow  conditions 
the  effects  of  turbulent  mixing  in  addition  to  compressibility  may  be  investigated. 

With  the  present  implementation  of  the  acetone  fluorescence  imaging  technique,  we  have 
managed  to  successfully  record  the  first  quantitative  images  of  the  separated  shear  layer. 
For  these  experiments,  the  separated  shear  layer  was  created  using  a  single-stream  shear 
layer  configuration.  The  results  obtained  from  the  single-stream  shear  layer  configuration 
confirm  the  expected  linear  growth  rate  in  the  far  field  as  seen  in  figure  5.  The  large 
size  of  the  pressure  vessel,  which  contains  the  test  section,  provides  ample  entrainment 
and  the  sharp  90-degree  corner  is  sufficiently  large  to  ensure  that  the  shear  layer  does 
not  reattach.  Images  similar  to  the  shear  layer  shown  in  figure  3  provide  quantitative 
information  on  the  instantaneous  spatial  structure  of  the  refractive  fluid  interfaces  in  this 
flow. 


Figure  5:  Preliminary  raw  image  of  the  separated  shear  layer  in  the  UCI  aero-optics 
Flow  Facility.  Flow  is  from  left  to  right,  and  the  high-speed  freestream  is  on  the  top  side. 
The  bottom  side  corresponds  to  the  ambient  air  in  the  pressure  vessel.  Linear  growth 
and  large-scale  features  are  evident  in  this  image  of  the  large-Reynolds-number  separated 
shear  layer.  The  dark  spots  correspond  to  residual  acetone  droplets  in  this  preliminary 
image.  These  have  been  eliminated  in  more  recent  experiments  in  our  facility  by  creating 
purely-gaseous  acetone  in  the  freestream  and/or  in  the  ambient  air  in  the  facility. 

The  present  flow-imaging  technique  is  currently  being  extended  to  record  directly  the 
density  field,  and  therefore  the  density  interfaces  generated  in  a  separated  shear  layer. 
This  procedure  extents  the  power  of  laser-induced  fluorescence  by  using  a  calibrated 
air/acetone- vapor  mixture  for  the  ambient  fluid  as  well  as  for  the  freestream  in  order 
to  capture  a  quantitative  representation  of  the  density  field.  Since  the  proposed  density 
whole-field  imaging  method  uses  purely  gaseous  fluid,  simultaneous  imaging  of  both  the 
flow  and  the  propagated  laser  wavefronts  is  enabled.  The  optical  wavefront  distortions 
will  be  conducted  using  a  high-resolution  spatial  Shack-Hartmann  sensor.  Initial  studies 
of  curved  separated  shear  layers  are  also  being  conducted  in  order  to  explore  the  use 
of  streamwise  flow  curvature  as  a  passive  means  for  flow  control  and  reduction  of  the 
aero-optical  distortions. 

The  facility  enables  the  simultaneous  imaging  of  the  density  field  and  the  optical  wave- 
fronts  through  five  10-inch  diameter  quartz  windows  which  allow  optical  access  to  the 
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interior  of  the  vessel.  A  new  imaging  technique  is  being  developed  where  laser-induced 
fluorescence  of  acetone  vapor,  seeded  in  the  flow  in  a  calibrated  way  to  remove  the  effects 
of  mixing  between  dissimilar  gas  streams,  can  be  used  to  provide  direct  measurements 
of  the  density  field.  The  concentration  of  acetone  vapor  is  regulated  and  introduced 
to  the  fully-enclosed  experiment  freestream  and  ambient  air  inside  the  test  section  us¬ 
ing  a  liquid  bubbling  technique.  The  acetone  vapor  concentration  was  maintained  in 
the  optically-thin  regime.  A  pulsed  Nd:YAG  laser  operated  at  a  wavelength  of  266  nm 
provides  the  excitation  source  for  acetone  which  yields  fluorescence  in  the  visible  range 
(350-550  nm).  The  first  flow  image  of  the  separated  shear  layer  recorded  in  the  new 
aero-optics  Flow  Facility  at  UC  Irvine  is  shown  in  figure  5.  One  of  the  main  objectives  of 
the  flow  imaging  technique  is  to  record  high  signal-to-noise  ratio  images  that  are  needed 
in  order  to  extract  information  of  the  refractive  interfaces  responsible  for  the  aero-optical 
distortions.  A  new  wavefront  sensor  is  being  developed  to  directly  measure  the  of  optical 
wavefronts  to  provide  high  spatial  resolution  of  the  optical  distortions  generated  by  the 
density  field.  It  is  a  custom  Hartmann  UV  sensor  based  on  a  DALSA  Piranha  4,096-pixel 
1,024-subaperture  detector.  This  will  permit  the  simultaneous  high-resolution  beam/flow 
imaging  using  the  new  aero-optics  flow  facility. 

4.  Conclusions  and  Implications 

The  Aero-Optics  Pressurized  Flow  Facility  at  UC  Irvine  is  a  unique  facility  that  is  es¬ 
pecially  designed  for  the  study  of  separated  flows.  It  permits  simultaneous  imaging  of 
both  the  flow  and  beam  behavior.  This  facility  is  expected  to  produce  archival  results  of 
aero-optical  interactions  and  to  provide  a  testing  environment  for  exploring  flow-control 
techniques  and  optical-control  techniques  to  reduce  aero-optical  distortions.  The  antici¬ 
pated  flow/beam  images  will  provide  a  high-resolution  database  which  can  be  analyzed 
to  determine  both  the  large-scale  and  small-scale  aero-optical  behavior  at  large  Reynolds 
numbers,  with  compressibility  effects,  using  the  new  Interfacial-Fluid-Thickness  frame¬ 
work  developed  by  the  PI  under  previous  AFOSR  support  (Catrakis  k,  Aguirre  2004 
[10],  AIAA  Journal,  accepted  for  publication;  Catrakis,  Aguirre,  &  Nathman  [11],  AIAA 
Journal ,  accepted  for  publication).1 
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Abstract 

A  new  approach  in  aero-optics  is  proposed  based  on  the  physical  thickness  of  refractive 
fluid  interfaces.  In  turbulent  flows,  particularly  at  large  Reynolds  numbers,  the  interfacial 
fluid  thickness  is  highly  variable.  The  role  of  this  interfacial  fluid  thickness  in  aero-optical 
interactions  is  examined  by  directly  relating  the  optical  path  length  (OPL)  to  the  interfacial- 
fluid-thickness  (IFT)  variations.  This  is  done  by  expressing  the  OPL  as  an  integral  of  the  IFT 
variations  along  the  beam  propagation  path.  We  demonstrate  the  proposed  IFT  approach 
on  refractive-field  measurements  in  large-Reynolds-number  (Re  ~  106)  high-compressibility 
(Mc  ~  1)  shear  layers  between  optically-different  gases.  Highly-irregular  networks  of  isolated 
high-gradient  interfaces  are  observed  at  various  transverse  locations  in  the  flow,  i.e.  both  in 
the  interior  and  near  the  outer  boundaries.  The  observation  that  the  high-gradient  interfaces 
are  spatially  isolated,  and  the  OPL  interpretation  in  terms  of  the  IFT  variations,  are  uti¬ 
lized  to  propose  and  demonstrate  a  new  modeling  approach  for  compressible  flow  where  the 
high-gradient  interfaces  are  the  dominant  elements  needed  to  capture  the  large-scale  optical 
distortions.  The  present  approach  suggests  a  new  point  of  view  for  relating  aero-optical 
distortions  to  the  flow  structure  in  terms  of  the  IFT  variations. 

1.  Introduction 

In  order  to  develop  techniques  for  the  prediction  and  control  of  aerooptical  phenomena 
in  turbulent  flows,  an  improved  understanding  is  needed  of  the  large-scale  properties  and, 
for  some  applications,  of  the  small-scale  properties  of  aerooptical  distortions  as  well  as  of  the 
turbulent  refractive  fluid  interfaces  that  generate  these  distortions  (e.g.  Jumper  &;  Fitzger¬ 
ald1  and  references  therein).  Practically,  the  large-scale  aerooptical  behavior  is  crucial  in  all 
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applications  involving  optical  beam  propagation  or  imaging  through  turbulent  shear  flows, 
such  as  the  flows  generated  by  airborne  vehicles.  The  small-scale  behavior  can  also  be  impor¬ 
tant  for  those  applications  that  require  high-resolution  and/or  long-range  optical  imaging  or 
beam  propagation.  In  addition,  it  is  important  to  know  in  practice  the  behavior  of  aeroop- 
tical  interactions  at  large  Reynolds  numbers.  Furthermore,  for  high-speed  flight,  it  is  also 
crucial  to  quantify  Mach-number  effects  on  the  aerooptical  behavior  (e.g.  Gilbert  &  Otten2 
and  references  therein). 

The  large-scale  organized  flow  behavior  is  known  to  provide  the  dominant  aerooptical  con¬ 
tributions  in  both  incompressible  turbulent  flows  (e.g.  Truman  &  Lee3;  Wissler  &  Roshko4; 
Dimotakis,  Catrakis,  &;  Fourguette5)  and  weakly-compressible  turbulent  flows  (e.g.  Jumper 
&  Fitzgerald1).  At  higher  flow  compressibilities,  however,  the  extent  to  which  the  large- 
scale  flow  behavior  is  organized  is  not  well  understood  (e.g.  Smits  &  Dussauge6;  Dimotakis, 
Catrakis,  &  Fourguette5).  The  large-scale  properties  of  aerooptical  distortions  at  high  com¬ 
pressibility,  and  the  manner  in  which  they  are  related  to  the  fluid-interfacial  behavior,  are 
also  not  well  understood. 

In  the  present  work,  we  propose  an  approach  to  examine  aerooptical  interactions  in 
terms  of  the  variations  in  the  thickness  of  refractive  fluid  interfaces.  Refractive  interfaces 
physically  exhibit  a  thickness  given  by  the  inverse  of  the  refractive-index  gradient  magnitude. 
This  interfacial-fluid-thickness  can  be  highly  nonuniform  in  turbulent  flows,  particularly  at 
large  Reynolds  numbers.  The  interfacial-fluid-thickness  variations  play  an  important  role 
in  aerooptics.  In  §2,  we  describe  the  interfacial-fluid-thickness  (IFT)  approach  and  show 
that  the  optical  path  length  (OPL)  can  be  expressed  directly  in  terms  of  the  variations  in 
interfacial  fluid  thickness.  In  §3,  we  demonstrate  this  approach  on  high-compressibility  large- 
Reynolds-number  shear  layers  of  optically-different  gases.  The  high-gradient  interfaces  are 
identified  and  their  behavior  is  compared  to  the  case  of  low  compressibility.  A  discussion  is 
included  on  the  issue  of  mixing  vs.  compressibility  effects.  In  §4,  we  propose  and  demonstrate 
a  new  modeling  approach  in  which  the  high-gradient  interfaces  are  used  to  produce  the  large- 
scale  aerooptical  distortions  at  high  compressibility.  In  a  companion  work,  we  examine  the 
small-scale  behavior  of  aerooptical  distortions  at  high  compressibility  and  large  Reynolds 
numbers,  by  means  of  a  new  method  for  characterizing  the  scale-local  structure  of  optical 
wavefronts  (Catrakis,  Aguirre,  &  Nathman7). 
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2.  Proposed  interfacial-fluid-thickness  (IFT)  approach  in  aerooptics 

Since  the  pioneering  aerooptics  studies  by  Liepmann8,9,  it  has  been  recognized  that  one 
of  the  central  goals  in  aerooptics  research  is  to  relate  the  optical-wavefront  distortions  to 
the  fluid-mechanical  behavior.  The  optical-wavefront  distortions  are  usefully  quantified  by 
the  optical  path  length  (OPL),  for  many  aerooptics  applications.  The  most  relevant  fluid- 
mechanical  quantity  in  aerooptics  is  the  refractive-index  field, 

n(x’t)  s  £  1  -  (1) 

where  c(x,  t)  is  the  local  speed  of  light  which  is  always  less  than  or  equal  to  the  speed 
of  light  Co  in  vacuum.  In  turbulent  flows,  the  refractive-index  field  n(x,  t)  can  be  highly 
nonuniform  particularly  at  large  Reynolds  numbers.  The  propagation  of  optical  wavefronts 
through  refractive-index  fields  is  governed  by  the  eikonal  equation,  which  can  be  expressed 
in  terms  of  the  OPL  as, 

|V  (OPL)|  =  n.  (2) 

The  eikonal  equation  is  useful  to  describe  fluid-optical  interactions  as  long  as  the  wavelength 
of  light  is  smaller  than  the  smallest  fluid-mechanical  (turbulent)  scale  (e.g.  Liepmann9)  and 
as  long  as  the  optical-beam  energies  are  low  enough  to  not  change  the  local  refractive  index. 

How  is  the  OPL  behavior  related  physically  to  the  structure  of  the  refractive-index  field 
and  in  particular  the  refractive  fluid  interfaces?  In  this  section,  we  present  a  framework  that 
is  useful  to  address  this  question  by  emphasizing  the  role  of  the  physical  thickness  of  the 
refractive  fluid  interfaces  in  determining  the  variations  in  the  OPL. 

As  a  starting  point,  we  recall  the  expression  used  for  the  definition  of  the  OPL  as  an 
integral  of  the  refractive  index  along  each  light  ray  in  geometric  optics  (e.g.  Dimotakis, 
Catrakis,  &  Fourguette5;  Jumper  &  Fitzgerald1,  and  references  therein),  i.e., 

OPL(x,f)  =  f  n(t,t)  df’,  (3) 

J  ray 

where  l  denotes  the  physical  distance  along  the  beam  propagation  path  for  each  light  ray. 
The  OPL  integral  in  equation  3  corresponds  to  inverting  the  eikonal  equation  2  for  the  OPL 
in  terms  of  the  refractive-index  field.  The  aerooptical  distortions  correspond  to  the  optical- 
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path  difference  (OPD)  given  by  OPD  (x,  t )  =  OPL  (x,  t)— OPLref  (x,  t ),  where  OPLref  (x,  t)  = 
/raynref(x,  t)d£  is  the  reference  OPL  that  would  correspond  to  the  undistorted  wavefronts, 
with  nref  denoting  a  reference  refractive  index,  e.g.  corresponding  to  freestream  conditions. 
The  optical  wavefronts  can  be  represented  as  isosurfaces  of  the  OPL,  i.e., 

OPL  (x,  t)  =  const. ,  (4) 

as  indicated  schematically  in  figure  1.  As  long  as  the  flow  speeds  are  small  relative  to  the 
speed  of  light,  and  the  optical-beam  propagation  distances  are  small  enough  for  light  to 
propagate  through  the  flow  before  it  has  evolved,  it  is  sufficient  to  think  of  the  OPL  integral 
in  equation  1  as  involving  only  the  instantaneous  spatial  structure  of  the  refractive-index 
field,  at  each  instant  in  time. 

As  the  optical  wavefronts  propagate  through  the  nonuniform  refractive-index  field  n(x,  t), 
the  aerooptical  interactions  physically  occur  across  the  refractive  fluid  interfaces.  These  are 
the  interfaces  on  which  the  refractive  index  n  is  constant.  It  is  important  to  understand 
the  role  of  these  interfaces.  While  these  fluid  interfaces  correspond  to  isosurfaces  of  the 
refractive-index  field,  i.e., 

n(x,  t )  =  const. ,  (5) 

it  is  crucial  to  recognize  that  the  refractive-fluid  interfaces  will  have  a  physical  thickness 
whereas  the  refractive-fluid  isosurfaces  are  geometrical  objects  with  zero  thickness.  It  is  the 
physical  thickness  of  the  interfaces  that  is  very  important  in  aerooptics.  We  can  introduce  the 
local  interfacial  fluid  thickness  hn,  defined  per  unit  n,  as  the  inverse  of  the  refractive-index 
gradient  magnitude,  i.e., 

'*"(x't)  s  m-  (6) 

The  distance  between  two  neighboring  isosurfaces,  corresponding  to  n  and  n  +  dn,  will  be 
hn  dn.  This  is  shown  schematically  in  figure  2,  for  a  single  interface,  and  in  figure  3  for  several 
neighboring  interfaces.  We  can  distinguish  between  the  general  case,  where  the  gradient 
magnitude  is  nonzero,  i.e.  |Vn|  >  0,  and  the  special  case  of  zero  gradient  magnitude, 
i.e.  |Vn|  =  0,  which  will  be  discussed  below.  In  regions  of  relatively-large  refractive- 
index  gradients,  the  isosurfaces  are  closely  spaced  and  the  interfaces  are  associated  with 
a  relatively-small  thickness.  In  regions  of  weak  refractive-index  gradients,  the  isosurfaces 
will  be  located  further  apart  and  the  interfaces  will  be  relatively  thicker.  The  interfacial 
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fluid  thickness  can  be  expected  to  be  highly  nonuniform  at  large  Reynolds  numbers.  This 
is  because  of  the  strongly-intermittent  character  of  fully-developed  turbulent  flows  which 
becomes  more  intermittent  with  increasing  Reynolds  number  for  both  incompressible  and 
compressible  flows  (e.g.  Sreenivasan10,  and  references  therein;  Smits  &  Dussauge6).  As 
discussed  by  Jumper  &  Fitzgerald1,  refractive-index  fluctuations  can  arise  in  pure  fluids, 
e.g.  where  density  fluctuations  are  induced  by  temperature  variations  in  low-speed  air  flows 
(e.g.  Jumper  &  Hugo11)  or  density  fluctuations  in  compressible  air  flows  (e.g.  Fitzgerald  & 
Jumper12),  or  in  mixtures  of  dissimilar  fluids  (e.g.  Brown  &  Roshko13;  Dimotakis,  Catrakis, 
&  Fourguette5) .  In  all  these  different  cases,  the  thickness  of  the  refractive  interfaces  can 
be  defined  in  the  same  manner,  i.e.  using  equation  6.  We  note  that  the  present  approach 
could  be  formulated  entirely  in  terms  of  the  refractive-index  gradient.  The  interfacial  fluid 
thickness,  however,  is  especially  helpful  for  physical  insight  into  the  aerooptical  interactions. 

Is  the  interfacial  fluid  thickness  finite?  Is  it  nonzero?  In  turbulent  flows,  the  interfa¬ 
cial  fluid  thickness  can  be  expected  to  be  both  nonzero  and  finite,  in  general,  on  physical 
grounds.  The  thickness  must  be  nonzero  wherever  the  local  refractive-index  gradient  mag¬ 
nitude  |Vn(x,  t)\  is  finite,  as  indicated  from  equation  6.  Only  an  infinite  gradient  can  lead 
to  a  zero  interfacial  fluid  thickness.  Physically,  it  is  clear  that  for  flows  of  real  fluids,  even  at 
large  but  finite  Reynolds  numbers,  the  finite  molecular  diffusivities  of  the  fluid  ensure  finite 
gradients  and  therefore  finite  interfacial  fluid  thicknesses.  A  related  observation,  also  indi¬ 
cated  from  equation  6,  is  that  the  thickness  must  be  finite  as  long  as  the  gradient  magnitude 
is  nonzero.  If  the  gradient  is  zero,  which  would  correspond  to  a  region  of  exactly-uniform 
refractive  index,  the  thickness  would  be  infinite  in  the  context  of  equation  6  and  one  then 
must  interpret  the  (infinity-times-zero)  product  hn  dri  as  a  Dirac  delta  function  whose  inte¬ 
gral  becomes  the  distance  given  by  the  extent  of  the  uniform-index  region  in  the  direction  of 
the  optical-ray  propagation.  In  summary,  we  can  expect  physically  that  the  interfacial  fluid 
thickness  must  be  finite  as  long  as  the  gradient  is  finite,  consistent  with  equation  6.  The 
thickness  is  an  interfacial  property,  therefore,  that  has  to  be  taken  into  account. 

Since  each  optical  ray  physically  propagates  through  refractive  interfaces,  one  can  intu¬ 
itively  expect  that  the  local  interfacial  fluid  thickness  should  determine,  at  least  in  part,  the 
local  contribution  to  the  OPL.  Can  this  be  seen  directly  in  the  OPL  integral  in  equation  3? 
This  can  be  done  by  rewriting  equation  3,  from  the  point  of  view  of  the  refractive  fluid 
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interfaces,  as 


OPL(x,t)  =  f  n(t,t)  hnye  \dn\ , 

J  ray 


(7) 


where  the  integration  is  now  performed  with  respect  to  the  refractive  index  n,  rather  than 
with  respect  to  the  spatial  distance  £,  and  hn/  is  the  effective  interfacial  fluid  thickness 
defined  as  the  component  of  the  interfacial  thickness  in  the  direction  of  optical  propagation, 
given  by 


hn,e  = 


|Vn|*  ’ 


(8) 


with  |  Vn|^  denoting  the  effective  gradient  magnitude,  i.e.  the  magnitude  of  the  component 
of  the  local  refractive-index  gradient  in  the  £  direction,  i.e.  in  the  direction  of  the  optical-ray 
propagation.  The  component  of  the  refractive-index  gradient  in  equation  7  is, 


|Vn|f 


\M 

d£  ’ 


(9) 


as  required  in  order  for  equations  2  and  6  to  be  consistent.  Since  the  refractive  index  n  could 
be  locally  increasing  or  decreasing  as  the  light  rays  propagate,  it  is  necessary  to  express  the 
differential  of  n  as  the  absolute- valued  differential  |dn|,  in  equation  7.  For  interfaces  locally 
normal  to  the  optical  rays,  the  gradient  component  |Vn|^  has  magnitude  identical  to  the 
magnitude  of  |Vn|.  Where  the  interfaces  are  locally  not  perpendicular  to  the  optical  rays, 
this  component  will  be  of  smaller  magnitude  than  |Vn|  and  the  effective  interfacial  fluid 
thickness  will  be  larger.  In  other  words,  the  effective  gradient  is  always  less  than  or  equal  to 
the  full  gradient,  i.e. 

|Vn|^  =  |Vn||cos0|  <  |Vn| ,  (10) 

and  the  effective  interfacial  fluid  thickness  is  always  greater  than  or  equal  to  the  full  interfacial 
fluid  thickness,  i.e. 

hn,e  =  j^-r  =  Msec0|  >  hn,  (11) 

|  vn|f 

where  the  angle  6,  taken  as  —  n  <  9  <  n,  quantifies  the  interfacial  orientation  relative  to 
the  optical-propagation  direction.  These  are  important  considerations  as  they  can  be  used 
to  develop  a  modeling  approach  based  on  the  interfacial  fluid  thickness,  described  below  in 
§4.  As  long  as  the  refractive-index  gradient  magnitude  is  not  zero,  i.e.  as  long  as  |Vn|  >  0, 
we  can  define  9  as  the  angle  between  the  refractive-index  gradient  vector  and  the  local 
optical-ray  propagation  vector.  The  refractive-index  gradient  vector  is  always  normal  to  the 
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local  refractive  interface.  Combining  equations  11  and  7,  and  writing  again  equation  7  for 
clarity,  we  see  that  the  OPL  integral  of  equation  3  can  be  expressed  directly  in  terms  of  the 
interfacial-fluid-thickness  variations  along  the  optical  beampropagation  path  as, 

OPL(x,  t)  =  f  n(£,t)  /i„|sec#|  |dn|  =  f  n{£,t)  hn^  |dn| ,  (12) 

J ray  J ray 

where  the  first  integral  is  in  terms  of  the  thickness  hn  and  relative  orientation  0  of  the  refrac¬ 
tive  interfaces,  and  the  second  integral  is  in  terms  of  the  effective  interfacial  fluid  thickness 
hnte.  As  mentioned  above,  in  the  context  of  the  interfacial  fluid  thickness  in  equation  5, 
these  integrals  require  that  |Vn|  ^  0.  In  those  regions  where  |Vn|  =  0,  i.e.  in  regions  of 
uniform  refractive  index,  the  (infinity-time-zero)  product  hndn  must  be  interpreted  again 
as  a  Dirac  delta  function  whose  integral  is  the  distance,  say  A£,  corresponding  to  the  extent 
of  the  uniform-index  region  in  the  optical-propagation  direction,  so  that  the  contribution 
to  the  OPL  integral  becomes  A  (OPL)  =  n  A£.  We  note  that,  in  uniform-index  regions, 
the  interfacial  orientation  6  has  no  meaning  and  is  not  needed.  In  those  instances  where 
the  fluid  interface  happens  to  be  locally  tangent  to  the  optical-propagation  direction,  i.e.  if 
9  =  ±  7t/2,  the  term  |  sec#|  will  be  infinite  but  in  such  cases  the  refractive  index  will  locally 
be  uniform  since  the  interface  will  be  aligned  with  the  optical-propagation  direction,  i.e. 
\dn\  =  0  in  such  cases.  In  those  cases,  therefore,  the  contribution  to  the  OPL  integral  will 
again  be  A  (OPL)  =  n  A£  with  A£  identified  as  the  length  of  the  interface  that  is  tangential 
to  the  optical-propagation  direction.  We  should  also  note  that,  in  general,  one  may  also  need 
to  take  into  account  other  possibilities  such  as  total  internal  reflection  or  the  development  of 
caustics.  These  possibilities  depend  on  the  magnitudes  of  refractive  gradients,  and  relative 
interfacial  orientations,  encountered  in  practice. 

In  summary,  the  proposed  interfacial-fluid-thickness  approach  is  based  on  relating  the 
OPL  to  the  interfacial-fluid-thickness  variations.  Whereas  the  integral  in  equation  3  is  con¬ 
ducted  over  space,  the  integrals  in  equations  7  and  12  are  expressed  as  integrals  over  the 
refractive  index  and  are  useful  to  determine  the  manner  in  which  the  refractive  interfaces 
physically  contribute  to  the  OPL.  In  addition  to  the  local  refractive  index  n,  equation  11 
shows  that  the  OPL  variations  arise  from  the  variability  in  the  interfacial  fluid  thickness 
hn  and  the  fluctuations  in  the  interfacial  orientation  9,  or  the  variations  in  the  effective 
interfacial  fluid  thickness  hn^.  Knowledge  of  the  variability  in  the  effective  interfacial  fluid 
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thickness,  and  its  relation  to  the  flow  dynamics,  can  be  expected  therefore  to  provide  physical 
insight  into  the  relation  between  the  OPL  behavior  and  the  interfacial  structure. 

3.  Demonstration  on  high- compressibility  turbulent-shear-layer  fluid  interfaces 

In  this  section,  we  demonstrate  the  use  of  the  interfacial-fluid-thickness  approach  on 
turbulent  high-compressibility  fluid  interfaces.  We  have  chosen  to  examine  the  applicabil¬ 
ity  of  the  proposed  approach  to  high-compressibility  fluid  interfaces  in  order  to  develop  a 
large-scale  aerooptics  modeling  methodology  useful  for  high-compressibility  flow  conditions. 
This  will  be  described  in  the  next  section  (§4).  There  are  available  large-scale  aerooptics 
modeling  methods  that  have  successfully  addressed  weakly-compressible  flows  (e.g.  Jumper 
&  Fitzgerald1)  and  incompressible  flows  (e.g.  Truman  &  Lee3;  Chew  &  Christiansen14; 
Dimotakis,  Catrakis,  &  Fourguette5).  These  methods  are  based  essentially  on  the  Brown- 
Roshko  large-scale  organized-structure  approach  (Brown  &  Roshko13)  with  direct  extensions 
to  weakly-compressible  flows,  and  will  be  discussed  briefly  in  §4.  It  is  known,  however, 
that  high-compressibility  turbulent  flows  exhibit  significant  differences  compared  to  weakly- 
compressible  or  incompressible  flows  (e.g.  Papamoschou  &  Roshko15;  Papamoschou  199116; 
Clemens  &  Mungal17).  A  more  general  approach  is  desirable,  to  be  able  to  address  aeroop- 
tical  effects  across  the  full  range  of  compressibilities,  and  it  is  toward  this  goal  that  the 
interfacial-fluid-thickness  approach  of  §2  can  be  particularly  useful. 

Figure  4a  shows  an  example  of  a  two-dimensional  spatial  streamwise  slice  of  the  refractive- 
index  field  in  a  shear  layer  between  optically-different  gases,  with  convective  Mach  number 
Mc  ~  1  and  Reynolds  number  Re  ~  106  based  on  the  large-scale  extent  Lg  of  the  flow 
(Dimotakis,  Catrakis,  &;  Fourguette5).  Also  shown  is  an  example  of  a  low-compressibility 
(Mc  ~  0.2)  shear-layer  image  in  figure  4b,  for  comparison.  These  flow  images  span  the  entire 
large-scale  transverse  extent  of  the  flow  and,  although  not  fully  resolved,  they  capture  a 
relatively-wide  range  of  scales  (~  500  :  1)  which  permits  a  study  of  the  large-scale  behavior. 
The  imaging  technique  relies  on  Rayleigh  scattering  and  on  using  dissimilar  gases  that  have 
a  large  difference  in  scattering  cross  section  and  a  relatively  large  difference  in  refractive 
index.  The  low-speed  freestream  contains  ethylene  in  both  the  high-  and  low-compressibility 
cases,  while  the  high-speed  freestream  contains  helium  or  nitrogen  for  the  high-  or  low- 
compressibility  cases  respectively.  The  freestream  refractive-index  difference  is  An  ~  6.6  x 
10~4  or  An  ~  4.0  x  10-4,  for  the  high-  or  low-compressibility  cases  respectively,  at  the  flow 
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conditions  investigated.  A  pulsed  laser  at  a  visible  wavelength  of  532  nm  and  with  energy  of 
~  300  mJ  per  pulse  provided  the  illumination  source  for  the  flow  imaging.  A  description  of 
the  flow  facility  and  further  details  of  the  experimental  technique  can  be  found  in  Dimotakis, 
Catrakis,  &  Fourguette.5 

It  should  be  noted  that  the  refractive-index  field  in  figure  4a  reflects  effects  of  both  mixing 
and  compressibility,  as  a  result  of  the  use  of  optically-different  gases,  and  may  not  directly 
represent  the  density-field  behavior,  as  has  been  discussed  by  Jumper  &  Fitzgerald.1  The 
issue  of  the  relative  effects  of  mixing  vs.  density  variations  is  important.  It  is  helpful  to  appre¬ 
ciate  first  that  nearly-uniform  concentration  regions,  in  the  case  of  mixing,  are  consistent  with 
and  directly  related  to  low-density  wells  corresponding  to  concentrated-vorticity  structures, 
e.g.  organized  vortical  structures.  In  fact  the  presence  of  nearly-uniform  regions  in  the  case 
of  mixing  provides  indirect  evidence  of  these  low-density  wells.  The  density-well  approach  of 
Jumper  &  Fitzgerald1  and  the  uniformly-mixed-region  approach  of  Dimotakis,  Catrakis,  & 
Fourguette5  are  both  related  to  the  original  Brown-Roshko13  organized-structures  approach 
since  the  Brown-Roshko  structures,  originally  recorded  at  incompressible  conditions,  are  the 
dominant  large-scale  vortical  structures  at  low  or  no  compressibility.  It  is  of  course  also 
important  to  remember  that  the  density  wells  become  deeper  as  the  level  of  compressibility 
rises.  Although  the  underlying  fluid  mechanics  is  the  same,  i.e.  the  same  vortical  struc¬ 
tures  generate  the  low-density  wells  and  the  uniformly-mixed  regions,  the  magnitude  of  the 
aerooptical  effects  associated  with  mixing  effects  vs.  pure  density  effects  can  be  expected 
to  be  quite  different.  It  is  clear  that  the  difference  between  these  two  effects,  aerooptically, 
can  be  large  depending  on  the  choice  of  the  dissimilar  gases  as  opposed  to  the  pure-air 
total-temperature-matched  freestreams.  Two  types  of  aerooptical  applications,  for  example, 
where  mixing  effects  can  be  important  are  high-speed-flow  applications,  where  a  film  coolant 
may  need  to  be  injected  near  the  optical  window  and  this  would  generate  mixing  between  the 
coolant  and  the  air,  i.e.  between  the  dissimilar  gases,  and  possible  future  aerooptics  appli¬ 
cations  whereby  one  may  need  to  be  able  to  vary  or  optimize  the  level  of  optical  distortions 
across  a  wide  range  of  magnitudes. 

The  primary  emphasis  of  the  present  work  is  on  the  interfacial-fluid-thickness  approach, 
which  is  a  general  approach  and  can  be  expected  to  be  useful  for  studying  both  mixing 
and  density  effects  in  the  future.  The  data  used  presently  to  demonstrate  this  approach 
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are  dominated  by  the  mixing  effect  but,  again,  the  main  approach  described  in  §2  is  of 
general  use,  i.e.  it  is  not  restricted  to  either  effect.  The  data  in  figure  4a  make  possible 
an  examination  of  the  interfacial-fluid-thickness  approach  at  high  levels  of  compressibility 
even  though  the  density-field  variations  associated  with  the  mixing  field,  can  be  expected 
to  exhibit  features  attributable  to  not  only  vortical  structures,  but  also  compression- wave 
regions  (e.g.  shocklets)  and  expansion-wave  regions.  The  refractive-index  field  in  figure  4a 
reflects  these  flow  mechanisms,  albeit  indirectly,  and  provides  the  opportunity  to  examine 
the  role  of  high-compressibility  fluid  interfaces  and  in  particular  their  physical  thickness  in 
generating  the  aerooptical  distortions. 

The  interfacial-fluid-thickness  fields,  hn,  corresponding  to  the  refractive-index  fields,  n 
of  figures  4a  and  b,  are  shown  in  figures  5a  and  b.  The  IFT  was  computed  by  first  evalu¬ 
ating  the  refractive-index  gradient  magnitude  |Vn|  and  subsequently  computing  its  inverse 
(Aguirre,  Ruiz-Plancarte,  &  Catrakis18),  according  to  the  IFT  definition  in  equation  6  of  §2. 
The  local  refractive-gradient  magnitudes,  and  interfacial  fluid  thicknesses,  are  necessarily 
under-/over-estimated  respectively  because  the  images  are  not  fully  resolved.  The  computed 
gradients  and  thicknesses  are  essentially  coarse-grained  values  at  the  image-resolution  scale. 
However,  these  values  permit  the  relative  spatial  variations  to  be  examined  and  this  is  what 
is  needed  in  the  present  context.  Since  the  refractive-index  data  of  figures  5a  and  b  are  from 
two-dimensional  slices  of  the  flow,  only  two  of  the  three  spatial  refractive-index  derivatives 
necessary  to  evaluate  the  full  gradient  magnitude  can  be  computed.  The  IFT  fields  in  fig¬ 
ures  5a  and  b  thus  correspond  to  the  in-plane  refractive-index  gradient  magnitude.  We  note, 
however,  that  only  the  component  of  the  full  gradient,  or  interfacial  fluid  thickness,  along  the 
optical-ray  propagation  direction  is  needed,  as  is  evident  in  equations  9,  10,  and  11  in  §2.  In 
other  words,  even  if  the  full  gradient  were  available,  one  would  compute  the  component  of  it 
in  the  optical-propagation  direction  and  this  would  be  equivalent  to  using  the  in-plane  com¬ 
ponent  in  addition  to  the  in-plane  interfacial  orientation,  for  optical  wavefronts  propagating 
in  the  plane  of  the  flow  shown  in  figures  5a  and  b,  or  figures  4a  and  b.  In  figures  5a  and  b, 
the  in-plane  gradient  magnitudes  are  shown  such  that  darker  regions  denote  higher  values 
of  the  gradient  magnitude.  Since  the  in-plane  interfacial  fluid  thickness  is  the  inverse  of  the 
in-plane  gradient  magnitude,  the  darker  regions  denote  locally- thin  high-gradient  interfaces. 
Also  evident  in  figures  5a  and  b  are  the  variations  in  the  interfacial  orientation.  Depending 
on  the  direction  of  optical  beam  propagation,  the  relative  interfacial  orientation  6  will  be 
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different  along  the  beam-propagation  path.  Figures  5a  and  b  indicate,  however,  that  the  in¬ 
terfacial  orientation  is  sufficiently  irregular  that,  at  least  locally,  there  is  a  wide  distribution 
of  interfacial  directions.  It  is  important  to  recall  the  meaning  of  equations  10  and  11  from 
§2,  i.e.  that  interfaces  of  a  given  thickness  will  have  an  effective  thickness  that  is  always 
greater  than  or  equal  to  it  for  any  optical-beam  propagation  direction.  In  other  words,  a 
given  gradient  magnitude  will  result  in  an  effective  gradient  magnitude  that  is  always  smaller 
than  or  equal  to  it,  irrespective  of  the  optical-beam  propagation  direction.  If  one  identifies, 
therefore,  the  regions  of  high-gradient  interfaces  in  the  flow,  then  there  can  be  no  effective 
gradient  magnitudes  that  are  higher  yet. 

Four  significant  observations  can  be  made,  at  high  compressibility,  on  the  basis  of  figure  5a 
and  other  interfacial-fluid-thickness  fields  at  the  same  flow  conditions: 

(a)  the  high-gradient  regions  occupy  a  relatively-small  part  of  the  turbulent  shear-flow  re¬ 

gion, 

(b)  the  high-gradient  regions  are  sheet-like  and  can  be  thought  of  as  high-gradient  interfaces, 

(c)  the  high-gradient  interfaces  exist  at  various  transverse  locations  in  the  shear-flow  region, 

and 

(d)  the  refractive-index  gradient  magnitude  exhibits  variations  along  the  high-gradient  in¬ 

terfaces. 

It  is  important  to  note  that  these  observations  refer  to  the  instantaneous  spatial  structure  of 
the  interfacial-fluid-thickness  field,  and  this  is  practically  what  is  most  relevant  for  aero  op¬ 
tics  since  it  is  the  instantaneous  flow  structure  that  needs  to  be  understood,  modeled,  and 
controlled.  The  ensemble-averaged  behavior  can  be  expected  of  course  to  be  quite  different 
and  simpler  with  smooth  ensemble-averaged  interfaces. 

Observation  (a)  refers  to  the  fact  that  the  high-gradient  regions  are  spatially  relatively 
isolated,  as  is  evident  in  figure  5a.  In  other  words,  large  parts  of  the  instantaneous  Bow 
region  have  relatively  low  refractive-gradient  magnitudes.  The  high-gradient  regions  appear 
to  occupy  only  a  small  fraction  of  the  flow  region.  This  is  consistent  with  computational 
results  at  low  Reynolds  numbers  in  studies  by  Samtaney,  Pullin,  &  Kosovic19  who  found 


11  of  31 


that  regions  of  high  gradient  magnitudes  in  isotropic  and  homogeneous  turbulence  occupy  a 
relatively-small  fraction  of  the  flow. 

Observation  (b)  indicates,  further,  that  the  high-gradient  regions  are  not  only  spatially 
isolated  but  axe  sheet-like,  i.e.  these  regions  are  confined  to  thin  layers  in  the  turbulent-flow 
region.  These  high-gradient  regions  can  be  thought  of  physically,  therefore,  as  locally-thin 
interfaces.  Sheet-like  structure,  in  three-dimensional  space,  is  consistent  with  the  string-like 
structure  evident  in  the  two-dimensional  spatial  images  or  slices  such  as  figure  5a. 

Observation  (c)  is  particularly  important  because  it  illustrates  a  major  qualitative  differ¬ 
ence  between  high-compressibility  flows  and  weakly-compressible  flows.  The  data  indicate 
that  the  high-gradient  interfaces,  although  spatially  isolated  and  occupying  a  relatively  small 
fraction  of  the  flow  region,  are  present  at  several  different  transverse  locations  in  the  instan¬ 
taneous  how  structure,  as  is  evident  in  the  high-compressibility  interfacial-fluid-thickness 
field  of  figure  5a  and  in  direct  contrast  to  the  low-compressibility  interfacial-fluid-thickness 
field  of  figure  5b.  Whereas  at  low  compressibility,  the  high-gradient  interfaces  are  mostly 
confined  to  the  instantaneous  outer  edges  of  the  flow,  it  is  clear  that  at  high  compressibility 
the  high-gradient  interfaces  can  be  found  both  in  the  interior  and  near  the  outer  bound¬ 
aries  in  instantaneous  realizations  of  the  flow.  This  is  in  sharp  contrast  with  the  behavior 
in  low-compressibility  turbulent  shear  flows  (Truman  &  Lee3;  Dimotakis,  Catrakis,  &  Four- 
guette5),  cf.  the  interfacial-fluid-thickness  behavior  in  figure  5b,  where  isolated  high-gradient 
regions  are  not  found  in  the  interior  of  the  flow,  as  a  result  of  large-scale  organized  structures 
which  confine  such  interfaces  mostly  to  the  outer  parts  of  the  shear  region.  At  the  high- 
compressibility  flow  conditions  examined  presently  in  figure  5a,  the  high-gradient  interfaces 
are  clearly  not  confined  to  the  outer  boundaries  of  the  flow  and  are  not  excluded  from  the 
interior  of  the  flow. 

Observation  (d)  is  also  important  because  it  indicates  that,  on  the  high-gradient  inter¬ 
faces,  the  refractive-index  gradient  is  not  constant,  i.e.  it  fluctuates  along  these  interfaces. 
This  must  be  taken  into  account  if  one  wishes  to  model  the  aerooptical  effects  of  these  in¬ 
terfaces,  and  this  will  be  addressed  in  the  proposed  modeling  approach  in  §4  below.  We 
note  that  the  variability  of  the  thickness  along  the  high-gradient  interfaces  is  attributable  to 
both  intermittency  effects  associated  with  the  large  Reynolds  number  and  the  presence  of 
compact  shocklets  associated  with  the  high  compressibility  of  the  flow. 
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These  observations,  and  the  dramatic  difference  between  the  high-compressibility  be¬ 
havior  and  the  low-compressibility  behavior,  are  illustrated  also  in  figures  6a  and  b.  The 
observation,  at  high  compressibility,  that  the  high-gradient  regions  are  spatially  isolated  even 
though  they  occupy  several  transverse  locations  in  the  flow,  and  the  interpretation  of  the 
optical  wavefront  phase  in  terms  of  the  interfacial  fluid  thickness  variations  (cf.  §2) ,  indicate 
that  these  isolated  regions  are  the  dominant  elements  for  the  large-scale  optical  distortions 
at  high  compressibility.  While  the  low-gradient  regions  are  of  wide  transverse  extent,  and 
therefore  do  contribute  to  the  OPL  integral,  it  is  the  structure  of  the  high-gradient  (locally- 
thin)  interfaces  that  can  form  the  basis  of  a  description  that  captures  the  large-scale  optical 
distortions,  as  explained  in  the  modeling  approach  proposed  in  the  following  section. 

4.  Proposed  large-scale  aerooptics  modeling  approach  based  on  high-gradient 
interfaces 

The  observations  of  the  instantaneous  structure  of  the  interfacial-fluid- thickness  field,  in 
§3,  and  the  interpretation  of  the  OPL  in  terms  of  the  interfacial-fluid-thickness  variations, 
in  §2,  can  be  used  to  develop  a  modeling  methodology  that  is  useful  for  capturing  the  large- 
scale  optical  distortions  at  high  compressibility.  Before  describing  the  proposed  modeling 
approach,  it  is  helpful  to  recall  the  previous  work  in  incompressible  and  weakly-compressible 
shear  layers  that  has  been  successful  in  modeling  the  large-scale  aerooptical  distortions. 
In  the  case  of  incompressible  shear  layers  of  dissimilar  gases,  the  Brown-Roshko  vortical 
structures  (Brown  &  Roshko13)  generate  large-scale  regions  of  nearly-uniform  refractive  index 
and  this  enables  the  use  of  the  outer  interfaces  to  model  the  large-scale  optical  distortions 
(Dimotakis,  Catrakis,  &  Fourguette5).  We  note,  however,  that  the  ’outer-interface’  model 
of  Dimotakis,  Catrakis,  &  Fourguette5  is  only  applicable  to  mixing  cases,  i.e.  it  does  not 
address  the  density  effects.  In  the  case  of  weakly-compressible  shear  layers,  there  are  large- 
scale  regions  of  reduced  density  and  pressure  in  the  flow  as  shown  by  Fitzgerald  &  Jumper.12 
These  large-scale  reduced-density  regions  are  essentially  weakly-compressible  analogues  of 
the  Brown-Roshko  structures.  Jumper  &  Fitzgerald1  have  developed  a  model  based  on  these 
large-scale  low-density  regions,  or  density  wells,  and  have  shown  that  it  is  able  to  reproduce 
well  the  large-scale  optical  distortions  at  low  compressibility.  The  present  IFT  framework 
(§2)  is  a  general  approach  and  thus  is  not  restricted  to  low  compressibilities. 

It  is  known,  from  various  investigations  of  high-compressibility  flows,  that  there  is  a 
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strong  qualitative  difference  between  the  structure  of  weakly-compressible  flows  and  high- 
compressibility  flows  (e.g.  Bogdanoff20;  Papamoschou  &  Roshko15;  Papamoschou16;  Samimy 
&  Elliott21;  Samimy,  Reeder,  &  Elliott22;  Clemens  &  Mungal17).  For  shear  layers,  in  addition 
to  the  well-known  reduction  in  the  growth  rate  there  is  a  substantial  difference  in  the  organi¬ 
zation  of  the  flow  and  in  particular  in  the  internal  structure  of  the  flow  as  the  compressibility 
level  is  increased.  In  particular,  accepting  the  convective  Mach  number  Mc  as  an  adequate 
measure  of  compressibility, 


Mc 


Ui  ~  U2 

Ol  +  0,2 


(13) 


where  U\^  and  a denote  the  freestream  speeds  and  sound  speeds,  respectively,  there  is 
evidence  that  the  flow  behavior  transitions  and  changes  significantly  at  convective  Mach 
numbers  above  a  certain  value  (Clemens  h  Mungal17).  Specifically,  in  the  range 


Me  >  0.6, 


(14) 


which  corresponds  to  high  compressibility  according  to  most  available  data,  the  turbu¬ 
lent  flow  structure  and  fluid  interfaces  appear  to  be  substantially  more  irregular,  and  in 
a  fundamentally-different  way,  compared  to  the  behavior  at  low  compressibility. 

Comparison  of  the  images  in  figures  5a  and  b,  example,  indicates  that  high-compressibility 
shear  layers  exhibit  multiple  regions  spanning  different  transverse  locations  where  the  refractive- 
index  gradients  are  relatively  weak.  At  the  boundaries  between  these  multiple  regions,  the 
high-gradient  interfaces  reside.  As  discussed  in  §3,  these  high-gradient  interfaces  appear  to 
form  a  highly-irregular  network  in  the  flow.  Is  there  a  way  to  model  the  interfacial-fluid- 
thickness  field  at  high  compressibility  in  order  to  reproduce  the  aerooptical  distortions?  In 
particular,  how  can  the  large-scale  optical  distortions  be  modeled  at  high  compressibility  and 
which  part  of  the  interfacial-fluid-thickness  structure  is  necessary  to  capture  the  large-scale 
optical  distortions? 

Our  proposed  modeling  approach,  at  high  compressibility,  is  based  on  the  crucial  observa¬ 
tion  in  §3  that  the  high-gradient  interfaces,  although  located  at  various  transverse  locations, 
are  spatially  isolated,  cf.  figures  5a  and  6a.  Because  the  high- gradient  interfaces  occupy  a 
relatively-small  fraction  of  the  total  turbulent-flow  region,  we  can  develop  a  modeling  ap¬ 
proach  with  these  interfaces  as  the  dominant  elements.  A  gradient-magnitude  threshold  is 
first  chosen  in  order  to  identify  the  high-gradient  interfaces.  We  note  that  it  is  crucial  to 
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retain  the  high-gradient  interfaces.  It  may  appear,  at  first  sight,  that  these  high-gradient 
interfaces  could  be  neglected  since  they  only  occupy  a  small  fraction  of  the  turbulent-flow 
region.  On  the  contrary,  the  high-gradient  interfaces  contribute  significantly  to  the  OPL 
integral  because,  even  though  the  thickness  hnj_  is  relatively  small,  the  differential  |dn|  is 
relatively  large  across  such  interfaces.  Significantly,  there  is  no  way  to  a  void  including  these 
high-gradient  interfaces  in  the  modeling  approach.  Fortunately,  however,  it  is  precisely  be¬ 
cause  these  interfaces  are  spatially  isolated  that  it  is  feasible  to  approximate  the  regions 
in  between  these  interfaces  as  zero-gradient  regions,  i.e.  as  uniform-index  regions.  This 
raises  the  important  question  of  how  to  propagate  the  beam,  in  the  model,  through  these 
regions.  For  comparison,  in  the  Dimotakis,  Catrakis,  &  Fourguette5  modeling  approach, 
the  refractive  index  in  the  large-scale  regions  bounded  by  the  outer  interfaces  is  modeled  on 
the  basis  of  the  value  predicted  from  the  large-scale  entrainment.  For  the  present  case  of 
high  compressibility,  however,  there  are  multiple  regions  through  which  the  beam  must  be 
propagated.  We  propose  that  the  way  to  do  this  is  to  use  the  gradient  information  along 
the  high-gradient  interfaces  to  update  the  refractive  index  as  the  OPL  integral  is  computed 
across  these  interfaces.  The  locations  of  the  high-gradient  interfaces  provide  the  length  scales 
needed  to  propagate  the  beam  through  the  regions  in  between  the  interfaces.  The  proposed 
modeling  approach  can  be  summarized  therefore  in  four  steps: 

(a)  the  high-gradient  interfaces  are  first  identified,  including  their  spatial  location  and  gra¬ 

dient  values, 

(b)  the  optical  beam  is  propagated  across  the  high-gradient  interfaces  by  updating  the  OPL 

integral  using  the  gradient  values  at  these  interfaces, 

(c)  the  regions  in  between  the  high-gradient  interfaces  are  modeled  as  zero-gradient  regions, 

and, 

(d)  the  value  of  the  refractive  index  in  each  zero-gradient  region,  between  the  high-gradient 

interfaces,  is  computed  using  the  gradient  value  on  the  interface  which  the  beam  prop¬ 
agates  across  as  it  enters  each  zero-gradient  region. 

A  schematic  illustrating  the  basic  idea  of  this  model  is  shown  in  figure  7  for  the  high- 
compressibility  case.  The  multiple  regions  in  the  interior  of  the  flow  are  denoted  as  B 
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through  G.  We  emphasize  that  it  is  important,  and  necessary,  to  retain  both  the  location  of 
and  gradient  values  along  the  high-gradient  interfaces  in  the  proposed  modeling  approach. 

The  utility  of  this  approach  can  be  demonstrated  with  the  high-compressibility  IFT  data 
examined  in  §3.  Figure  8  shows  a  modeled  interfacial-fluid-thickness  field  where  ~  50%  of 
the  refractive-gradient  magnitudes  have  been  set  to  zero,  compared  to  the  original  field  in 
figure  5a.  The  OPL  integral  was  computed  using  this  modeled  IFT  field  by  following  the 
above  four  steps.  The  comparison  between  the  full  aerooptical-wavefront  OPL  (solid  curve) 
and  the  modeled  OPL  (dashed  curve)  in  shown  in  figure  9.  Good  agreement  is  evident 
in  terms  of  the  large-scale  OPL  variations.  Different  thresholds  were  examined  spanning 
almost  the  entire  range,  i.e.  from  almost  100%  to  almost  0%.  As  expected,  the  ability  of 
the  thresholded  field  to  capture  the  large-scale  wavefront  distortions  was  found  to  decrease 
with  decreasing  threshold.  The  important  finding  is  that  with  modest  threshold  levels  one 
can  capture  the  large-scale  wavefront  signature.  The  50%  threshold  employed  in  figure  8 
was  chosen  as  an  example  to  illustrate  that  the  high-gradient  interfaces  contain  the  essential 
information  to  model  the  large-scale  optical  distortions. 

Figure  9  indicates  that  the  proposed  modeling  approach  can  be  used  to  reproduce  the 
large-scale  optical  distortions  at  high  compressibility,  in  a  shear  layer,  using  the  high-gradient 
information.  High-gradient  interfaces  both  in  the  interior  and  near  the  outer  boundaries 
of  the  flow  are  essential  to  capture  the  large-scale  wavefront  distortion.  However,  instead 
of  the  full  interfacial-fluid-thickness  field,  information  on  the  high-gradient  interfaces  only 
appears  to  be  sufficient  to  model  the  large-scale  aerooptical  distortions.  The  location  of  the 
high-gradient  interfaces  and  the  value  of  the  gradient  across  (or  of  the  thickness  of)  these 
interfaces  is  enough  to  capture  the  dominant  contributions  that  generate  the  large-scale 
aerooptical  distortions.  The  optical  wavefronts  are  propagated  through  the  high-gradient 
(locally-thin)  interfaces  and  the  gradient  value  across  these  interfaces  is  used  to  compute 
the  difference  in  the  optical  path  length,  while  in  the  low-gradient  regions  between  the 
high-gradient  interfaces  the  wavefronts  are  propagated  neglecting  the  presence  of  the  low 
gradients,  i.e.  as  if  those  regions  are  zero-gradient  regions.  The  high-gradient  interfaces,  in 
this  example,  are  identified  in  an  absolute  sense,  i.e.  by  thresholding  the  entire  instantaneous 
gradient  field.  A  further  refinement,  potentially,  could  be  envisaged  by  retaining  those 
interfaces  that  have  high  gradients  relative  to  their  neighboring  interfaces.  The  proposed 
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modeling  approach  offers  a  reduction  in  the  amount  of  flow  information  needed  to  capture 
the  large-scale  aerooptical  distortions  at  high  compressibility. 

5.  Conclusions  and  implications 

A  new  interfacial-fluid-thickness  (IFT)  approach  is  proposed  to  examine  optical-wavefront 
propagation  through  turbulent  flows  in  terms  of  the  variability  in  the  physical  thickness  of 
the  refractive  fluid  interfaces.  The  interfacial  fluid  thickness  is  given  by  the  inverse  of  the 
refractive-index  gradient  magnitude,  and  it  is  highly  variable  at  large  Reynolds  numbers 
and  high  compressibility.  The  IFT  plays  an  important  role  in  aerooptics  as  indicated  by 
expressing  directly  the  optical  path  length  (OPL)  in  terms  of  the  interfacial-fluid-thickness 
variations.  Demonstration  of  this  approach  on  refractive-field  data  in  shear  layers  at  high 
compressibility  (Mc  ~  1)  and  large  Reynolds  number  (Re  ~  106)  reveals  that  regions  of 
high  refractive-gradient  magnitudes  are  located  at  various  transverse  locations  in  the  flow, 
i.e.  both  in  the  interior  and  near  the  outer  boundaries,  and  form  highly-irregular  networks. 
Significantly,  the  high-gradient  regions  are  found  to  be  spatially  isolated,  i.e.  occupy  a 
relatively-small  fraction  of  the  shear  layer.  This  observation,  coupled  with  the  interpretation 
of  the  OPL  in  terms  of  interfacial-fluid-thickness  variations,  are  utilized  to  propose  and 
demonstrate  a  new  modeling  approach  where  the  high-gradient  interfaces  are  the  dominant 
elements  necessary  to  reproduce  the  large-scale  optical  distortions  at  high  compressibility. 
Both  the  location  of  and  the  gradient  values  along  the  high-gradient  interfaces  are  utilized  in 
this  modeling  approach.  As  an  example,  a  reduction  of  ~  50%  in  the  amount  of  interfacial 
information  needed  is  shown  to  reproduce  well  the  large-scale  optical  distortions  at  high 
compressibility.  The  present  results  suggest  applicability  of  the  proposed  approach  to  other 
high-compressibility  flows,  by  modeling  the  large-scale  optical  distortions  in  terms  of  the 
structure  of  the  high-gradient  interfaces.  The  proposed  approach  can  be  expected  to  be 
useful  in  studies  of  the  dynamics  of  aerooptical  interactions  at  high  compressibility  (e.g. 
Thurow,  Samimy,  Lempert,  Harris,  Widiker,  k  Duncan23),  in  computational  modeling  and 
simulations  of  aerooptical  interactions  in  high-speed  flows  (e.g.  Jones  k  Bender24),  as  well 
as  in  efforts  to  extend  fluid-mechanical  techniques  that  control  or  regularize  aerooptical 
interactions  in  weakly-compressible  flows  (e.g.  Stanek,  Sinha,  Seiner,  Pearce,  k  Jones25), 
Siegenthaler,  Jumper,  k  Asghar26)  to  higher  flow  compressibilities. 
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Figure  1.  Schematic  of  a  refractive-index  field  n(x,y;z,t)  at  high  compressibility,  shown  as 
the  gray-level  image  at  the  inset,  with  flow  from  left  to  right.  Also  shown  is  a  planar  incident 
optical  wavefront  (top),  propagating  from  top  to  bottom,  and  a  propagated  distorted  optical 
wavefront(bottom).  Ls  is  the  large-scale  transverse  extent  of  the  flow  and  measures  the 
extent  of  the  outermost  edges  of  the  turbulent  refractive  fluid  interfaces. 

Figure  2.  Schematic  representation  of  a  refractive  interface,  its  interfacial  fluid  thickness 
/in(x,  f),  and  two  neighboring  isosurfaces  of  the  refractive-index  field  n(x,  t). 

Figure  3.  Schematic  of  several  refractive  fluid  interfaces  and  the  variability  in  interfacial 
fluid  thickness. 

Figures  4a  and  b.  (a)  Refractive-index  field  in  a  streamwise  slice  of  a  high-compressibility 
(Mc  ~  1)  large-Reynolds-number  (Re  ~  106)  shear  layer  between  optically-different  gases.5 
(b)  Low-compressibility  (Mc  ~  0.2)  shear-layer  refractive  field.5 

Figures  5a  and  b.  (a)  Interfacial-fluid-thickness  field  hn  at  high  compressibility  corresponding 
to  the  refractive-index  field  of  figure  4.  Darker  regions  denote  thinner  interfaces,  i.e.  the 
dark  regions  correspond  to  the  high-gradient  interfaces,  (b)  Low-compressibility  interfacial- 
fluid-thickness  field. 

Figure  6a  and  b.  (a)  Schematic  of  network  of  spatially-isolated  high-gradient  interfaces  at 
high  compressibility.  Darker  regions  denote  higher  gradients  or  locally-thinner  interfaces, 
(b)  Schematic  of  high-gradient  interfaces  at  low  compressibility. 

Figure  7.  The  labeled  regions  in  the  interior  of  the  flow,  i.e.  regions  B  through  G,  correspond 
to  the  low-gradient  regions  at  high  compressibility.  In  the  proposed  modeling  approach, 
the  high-gradient  interfaces,  or  boundaries  between  the  low-gradient  regions,  are  treated 
as  the  dominant  elements  necessary  to  model  the  large-scale  optical  distortions  at  high 
compressibility. 

Figure  8.  Modeled  interfacial-fluid-thickness  field,  with  ~  50%  of  the  gradient  magnitudes 
neglected,  at  high  compressibility,  cf.  figure  5a. 

Figure  9.  Comparison  of  the  wavefront  OPD,  at  high  compressibility,  between  the  full 
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wavefront  OPD  (solid  curve)  using  the  complete  interfacial-fluid- thickness  field  (figure  5a), 
and  the  modeled  OPD  (dashed  curve)  using  the  high-gradient  interfaces  with  ~  50%  of  the 
gradient  magnitudes  neglected  (figure  8).  The  spatial  and  OPD  coordinates  are  normalized 
based  on  the  large-scale  flow  extent  L$,  in  a  similar  manner  as  in  Dimotakis,  Catrakis,  k. 
Fourguette5.  The  large-scale  optical  distortions,  at  high  compressibility,  can  be  captured 
in  terms  of  the  high-gradient  interfaces  only,  even  with  a  ~  50%  reduction  in  interfacial 
information,  demonstrating  the  utility  of  the  proposed  IFT  modeling  approach. 
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Abstract 

A  new  technique  is  proposed  which  enables  the  scale-local  characterization  of  aerooptical 
wavefronts.  Because  optical  wavefronts  degraded  by  turbulent  flows  are  physically  highly 
anisotropic  and  exhibit  distortions  over  a  wide  range  of  scales,  a  method  is  needed  to  ex¬ 
amine  the  wavefront  structure  at  varying  degrees  of  anisotropy  and  as  a  function  of  scale. 
We  introduce  an  aerooptical-wavefront-anisotropy  (AWA)  parameter  as  the  ratio  of  scaling 
factors  for  the  wavefront  distortions  and  for  the  wavefront  transverse  extent.  This  AWA 
parameter,  combined  with  box  counting,  enables  a  scale-local  anisotropic  examination  of 
the  wavefronts.  We  demonstrate  this  technique  on  wavefronts  derived  from  experiments  on 
large-Reynolds-number  (Re  ~  106)  high-compressibility  (Mc  ~  1)  turbulent  shear  layers 
between  optically-different  gases.  Variation  of  the  AWA  parameter  and  scale-local  exam¬ 
ination  of  the  distortions  reveal  the  presence  of  anisotropic  self-similarity,  or  self-affinity, 
at  small  scales  spanning  nearly  a  decade.  This  finding  shows  the  utility  of  the  technique 
to  detect  scaling  in  large-Reynolds-number  flow  experiments  from  the  aerooptical  behavior. 
The  present  finding  and  technique  provide  the  key  ingredients  to  extrapolate  the  small-scale 
properties  of  compressible-flow  aerooptical  wavefronts  to  higher  Reynolds  numbers,  and  are 
also  useful  for  modeling  and  for  computational  aerooptics. 

1.  Introduction 

The  large-scale  properties  of  aerooptical  interactions  in  turbulent  flows  can  be  expected 
to  be  strongly  dependent  on  the  flow  geometry  and  boundary  conditions  (e.g.  Truman  & 
Lee1;  Dimotakis,  Catrakis,  &  Fourguette2;  Jumper  &;  Fitzgerald3  and  references  therein).  In 
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contrast,  at  least  at  large  Reynolds  numbers  and  for  incompressible  flow  conditions,  there  is 
evidence  that  the  small-scale  flow  structure  can  exhibit  scaling  or  self-similar  behavior  which 
is  weakly,  if  at  all,  sensitive  to  the  flow  geometry,  in  the  sense  that  small-scale  turbulence 
statistics  appear  to  be  independent  of  the  flow  geometry  or  large-scale  behavior  to  a  good 
approximation  (e.g.  Sreenivasan4;  Catrakis,  Aguirre,  Sc  Ruiz-Plancarte5).  The  Reynolds 
number  needs  to  be  sufficiently  large  for  the  turbulent  behavior  to  be  fully-developed.  This 
requires  at  a  minimum  that  the  Reynolds  number  Re  based  on  the  large  scales  must  be  at 
least  Re  >  104  (e.g.  Dimotakis6).  In  airborne  applications  requiring  high-resolution  long- 
range  optical  imaging  or  beam  propagation,  which  are  necessarily  associated  with  large 
Reynolds  numbers,  a  complete  description  requires  knowledge  of  both  the  small-scale  prop¬ 
erties  and  the  large-scale  behavior. 

A  useful  point  of  view  in  developing  techniques  to  characterize  as  a  function  of  scale 
highly-irregular  objects,  such  as  optical  wavefronts  or  fluid  interfaces  in  turbulence,  is  to  ex¬ 
amine  geometric  properties  and  particularly  the  extent  to  which  self-similar/fractal  behavior 
(e.g.  Sreenivasan4;  Villermaux  Sc  Innocenti7)  or  scale-dependent  behavior  (e.g.  Catrakis  Sc 
Dimotakis8)  is  exhibited  at  the  small  scales.  Since  the  structure  of  aerooptical  wavefronts  is 
essentially  determined  by  the  turbulent  flow  structure,  one  may  expect  that  the  aerooptical 
distortions  will  exhibit  self-similar  behavior  at  small  scales  if  the  refractive  fluid  interfaces  are 
self-similar.  This  has  been  suggested  theoretically  for  optical  wavefronts  distorted  by  incom¬ 
pressible  turbulent  flows  (e.g.  Hentschel  Sc  Procaccia9;  Schwartz,  Baum,  Sc  Ribak10,11).  For 
compressible  flows,  there  is  some  evidence  that  the  small-scale  structure  of  fluid  interfaces 
can  exhibit  self-similarity  (e.g.  unpublished  work  of  Sreenivasan  Sc  Johnson  cited  in  Smits 
Sc  Dussauge12;  Poggie13).  Such  self-similarity  can  expected  to  be  at  most  weakly  dependent 
on  the  large-scale  flow  properties.  Thus,  for  compressible  turbulence,  one  might  expect  self¬ 
similar  scaling  behavior  for  the  small-scale  properties  of  both  the  turbulent  fluid  interfaces 
and  the  turbulence-distorted  optical  wavefronts.  Self-similar  behavior  of  the  optical  wave- 
fronts,  if  found  at  high  compressibility,  would  permit  extrapolating  the  small-scale  optical 
behavior  to  the  large  Reynolds  numbers  of  interest  in  practical  applications  in  high-speed 
flight. 

In  the  present  work,  we  introduce  an  aerooptical-wavefront-anisotropy  (AWA)  parame¬ 
ter  defined  as  a  ratio  of  scaling  factors  for  the  optical  path  difference  (OPD)  and  for  the 
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wavefront  spatial  extent.  This  is  described  in  §2  in  the  context  of  the  physical  anisotropy 
of  aerooptical  wavefronts.  The  AWA  parameter  enables  the  examination  of  the  small-scale 
geometrical  properties  of  the  wavefronts  using  an  anisotropic  extension  of  the  box  count¬ 
ing  method  and  this  is  developed  in  §3.  We  demonstrate  the  utility  of  this  technique 
on  wavefront  data  derived  from  large-Reynolds-number  high-compressibility  shear  layers 
between  optically-different  gases.  A  description  of  the  experimental  technique,  and  a  dis¬ 
cussion  of  the  important  issue  of  mixing  vs.  density  effects  in  aerooptics,  can  be  found 
in  the  companion  paper  (Catrakis  &;  Aguirre14)  and  references  therein.  In  §4  the  results 
of  the  scale-local  examination  of  the  wavefront  data  are  presented.  The  AWA  parameter 
is  found  to  be  useful  to  discern  anisotropically  self-similar,  or  self-affine,  behavior  of  the 
OPD  at  small  scales.  This  scaling  behavior  can  be  expected  in  other  large-Reynolds-number 
highly-compressible  flows  since  the  small-scale  behavior  is  at  most  weakly  dependent  on  the 
large-scale  behavior,  as  long  as  the  Reynolds  number  is  large  enough  for  the  turbulent  flow 
to  be  fully-developed.  Consequences  of  the  present  findings  are  discussed,  in  the  conclu¬ 
sions,  regarding  large-Reynolds-number  experimental  studies,  modeling,  and  computational 
aerooptics.  Considerations  of  the  large-scale  aerooptical  distortions,  using  a  new  interfacial- 
fluid- thickness  approach,  are  the  subject  of  the  companion  paper  (Catrakis  &  Aguirre14). 

2.  Aerooptical- wavefront-anisotropy  (AWA)  parameter  in  large-Reynolds-number 
turbulent  flows 

One  of  the  practical  as  well  as  fundamental  objectives  in  aerooptics  is  to  develop  the 
means  to  extrapolate  the  behavior  of  the  optical-wavefront  structure  to  larger  Reynolds 
numbers,  including  effects  of  compressibility.  Knowledge  of  the  small-scale  structure  of 
aerooptical  wavefronts  is  needed  toward  this  objective.  Optical  wavefronts  propagating 
through  turbulent  flows  can  be  expected  to  exhibit  behavior  that  is  similar  to,  or  related 
closely  to,  the  flow  structure.  To  be  able  to  extrapolate  the  behavior  of  a  turbulence¬ 
generated  refractive-index  field  n(x,  t)  to  large  Reynolds  numbers,  it  is  important  that  the 
Reynolds  number  of  the  flow  under  study  be  at  least  above  the  minimum  value  associated 
with  the  transition  to  fully-developed  turbulence  (e.g.  Dimotakis6),  i.e. 

Re  >  104.  (1) 

As  long  as  the  turbulent  flow  is  fully  developed,  i.e.  satisfies  at  least  the  criterion  in  equa- 
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tion  1,  the  range  of  scales  exhibited  by  the  refractive  field  n(x,t )  increases  with  increasing 
Reynolds  number,  as  is  well  known,  as, 

^  ~  Re3/4,  (2) 

^min 

for  gas-phase  flows  (Kolmogorov15),  where  Amjn  is  the  smallest  scale  or  Kolmogorov  scale 
and  Amax  is  the  largest  scale  identified  as  the  extent  of  the  flow  region  through  which  the 
wavefronts  propagate  or  alternatively  as  the  extent  of  large-scale  shear  in  the  flow.  An 
important  issue,  in  the  context  of  the  present  work,  is  that  there  is  a  narrower  range  of 
scales,  which  appears  only  for  large-Reynolds-number  flows  and  refractive-index  fields,  and 
this  is  the  range  of  scales  bounded  from  above  by  the  Liepmann-Taylor  scale  Alt  and  from 
below  by  Amin,  i.e., 

Amin  <  A  <  Alt  ~  Amax  Re-1/2 ,  (3) 

where  similarity  properties  (e.g.  Kolmogorov15),  including  possibly  self-similarity  (e.g.  Sreeni- 
vasan4),  are  present.  The  Liepmann-Taylor  scale  Alt  denotes  the  largest  of  the  small  scales 
that  are  not  directly  affected  by  large-scale  effects  or,  equivalently,  the  smallest  scale  that  is 
directly  affected  by  large-scale  effects.  Such  properties  are  important  because  they  can  be 
expected  to  help  extrapolate  aerooptical  behavior  to  larger  Reynolds  numbers,  encountered 
in  practice,  in  both  incompressible  and  compressible  flows  (e.g.  Smits  k  Dussauge12). 

In  order  to  characterize  aerooptical  distortions  over  the  wide  range  of  scales  associated 
with  large  values  of  the  Reynolds  number  (cf.  equations  1-3),  techniques  are  needed  that 
enable  the  description  of  optical- wavefront  structure  as  a  function  of  scale.  Most  previous 
studies  of  the  small-scale  flow  behavior  and  of  the  optical- wavefront  behavior  have  focused  on 
spectral  descriptions  or,  equivalently,  on  descriptions  based  on  second-order  structure  func¬ 
tions  (e.g.  Andrews  k  Phillips16;  Dimotakis,  Catrakis,  k,  Fourguette2).  Spectral  descriptions, 
however,  characterize  the  aerooptical  behavior  in  Fourier  space  and  do  not  uniquely  iden¬ 
tify  the  geometrical  structure  of  the  turbulent  flow  or  optical  wavefronts.  This  is  because 
power  spectra  retain  no  phase  information.  The  description  of  the  geometrical  structure 
of  turbulent  flows  in  physical  space  is  a  relatively-recent  development  (e.g.  Sreenivasan4; 
Catrakis,  Aguirre,  k  Ruiz-Plancarte  5)  and  relies  on  box-counting  techniques.  A  generalized 
version  of  the  box-counting  technique,  that  can  characterize  anisotropic  surfaces  such  as 
optical  wavefronts,  will  be  proposed  and  described  in  §3.  In  the  present  section,  the  issue 
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of  optical-wavefront  anisotropy  will  be  first  addressed  as  it  is  a  crucial  physical  aspect  of 
wavefront  structure  and  an  important  ingredient  of  the  new  box-counting  technique  of  §3. 

The  governing  equation  for  the  propagation  of  optical  wavefronts,  at  the  level  of  de¬ 
scription  of  geometrical  optics,  is  the  eikonal  equation  for  the  optical  path  length  (OPL)  or, 
equivalently,  the  inverse  or  solution  to  the  eikonal  equation  which  is  the  integral  for  the  OPL 
in  terms  of  the  refractive  index  along  the  optical  beam  propagation  path,  i.e. 

| V  (OPL) |  =  n,  or,  OPL(x,t)  =  f  n(£,t)  d£  =  f  n(£, t)  hn^  |dn| ,  (4) 

J  ray  J  ray 

where  the  first  integral  is  over  the  physical  distance  i  along  the  propagation  path.  In  the 
second  integral,  the  OPL  is  expressed  directly  in  terms  of  the  effective  refractive  interfacial 
fluid  thickness  (IFT)  which  is  hr,x  =  l/|Vn|f  (Catrakis  &  Aguirre14).  As  usual,  the  refractive 
index  is  n  =  co/c,  where  Co  is  the  speed  of  light  in  vacuum  and  c  is  the  local  speed  of  light. 
Since  the  velocity  of  light  is  related  to  the  propagation  distance  through  c  =  d£/dt,  we  can 
also  express  the  OPL  integral  as  OPL  =  J  nd£  =  cq  J  (l/c)d£  =  Co  J dt,  i.e.  it  follows 
that  the  OPL  integral  is  essentially  the  integrated  propagation  time,  or  wavefront  travel 
time,  even  though  the  integral  is  originally  defined  as  a  spatial  integral. 

In  the  present  work,  the  emphasis  is  on  the  optical-wavefront  structure.  The  optical 
wavefronts  are  physically  the  isosurfaces  of  the  OPL,  i.e. 

OPL  (x,  t )  =  const. .  (5) 

This  simple  and  important  statement  is  clear  physically  because,  as  stated  above,  the  OPL 
integral  is  essentially  a  propagation-time  integral,  i.e.  (OPL)/co  =  f  dt.  A  surface  of 
constant  OPL  is,  therefore,  a  surface  on  which  all  points  have  travelled  over  the  same  propa¬ 
gation  time,  ie.  duration,  and  therefore  this  surface  is  the  wavefront.  By  simply  normalizing 
the  OPL,  one  can  see  that  the  OPL  integral  is  also  essentially  an  optical-phase  integral.  We 
see,  therefore,  that  isosurfaces  of  the  OPL,  or  of  the  travel  time,  or  of  the  optical  phase, 
are  all  equivalent  physical  representations  of  the  optical  wavefronts.  It  is  noteworthy  to 
recall,  in  contrast,  that  the  wavefront-location  integral  /  cdt  =  J  d£  does  not  of  course  have 
this  property,  i.e.  its  isosurfaces  do  not  represent  the  wavefronts,  since  the  isosurfaces  of 
the  wavefront-location  integral  are  just  surfaces  of  constant  propagation  distance.  This  il- 
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lustrates  what  is  perhaps  the  most  important  physical  property  of  the  OPL  integral,  i.e. 
that  the  OPL  integral  is  physically  and  intuitively  very  helpful  precisely  because  it  has  the 
property  that  the  OPL  isosurfaces  are  the  wavefronts  themselves. 

In  turbulent  flows,  and  especially  at  large  Reynolds  numbers,  the  OPL  isosurfaces  or 
wavefronts  can  be  expected  to  be  highly  irregular.  In  the  present  work,  our  objective  is  to 
characterize  the  physical  structure  of  the  optical  wavefronts  as  a  function  of  scale  allowing 
for  anisotropic  scaling  behavior.  The  first  issue  that  needs  to  be  addressed  is  that  aerooptical 
wavefronts  are  physically  highly-anisotropic  surfaces.  This  is  because,  without  any  normal¬ 
ization,  the  spatial  transverse  extents  of  the  optical  wavefronts  are  in  practice  much  larger 
than  the  OPL  variations.  This  is  especially  so  because  the  refractive-index  variations,  e.g.  as 
quantified  by  the  root-mean-square  magnitude,  are  relatively  small  when  compared  to  unity 
for  gas-phase  flows.  Consequently,  as  the  integrals  in  equation  4  suggest,  the  OPL  variations 
will  be  much  smaller  than  the  spatial  transverse  extent  (i.e.  aperture)  of  any  optical  beam  in 
practice.  The  OPL  variations  are  best  expressed  through  the  optical  path  difference  (OPD) 
where  OPD  =  OPL  —  fray  nref(x,  t)d£,  where  nref  denotes  a  reference  refractive  index,  e.g. 
corresponding  to  a  freestream  condition.  Physically,  therefore,  the  wavefront  anisotropy  is 
present  because, 

|An|rms  <1,  and,  OPDrms  <  \A£\ ,  (6) 

where  |An|rms  and  OPDrms  denote  respectively  the  root-mean-square  magnitudes  of  the 
refractive-index  variations  and  of  the  OPD  of  the  distorted  wavefront,  and  | At |  denotes  a 
measure  of  the  total  optical-beam-propagation  path  which  is  |A^|  ~  L  or  Amax,  i.e.  com¬ 
parable  to  the  large-scale  extent  L  of  the  flow  as  far  as  the  aerooptical  interactions  in  the 
turbulent  region  are  concerned.  This  strong  anisotropy  of  aerooptical  wavefronts  refers  to 
their  physical,  unnormalized  structure.  One  may  normalize  the  OPD  so  that  the  effective 
degree  of  anisotropy  can  change.  By  also  scaling  the  spatial  coordinates,  one  can  examine 
the  wavefronts  at  varying  degrees  of  anisotropy,  as  will  be  shown  below  and  in  §3. 

A  schematic  showing  a  planar  wavefront  incident  on  a  three-dimensional  optically-active 
turbulent-flow  region  and  the  irregular  wavefront  emerging  from  the  flow  in  three  dimensions 
is  depicted  in  figure  1.  The  distorted  wavefront  is  highly  irregular  and  highly  anisotropic.  We 
can  also  consider  the  OPL  or  OPD  profiles,  which  are  themselves  highly  irregular.  The  ideas 
below  are  applicable  to  both  the  wavefront  surfaces  and  the  OPL  or  OPD  profiles.  Taking 
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the  x-  and  y-directions  as  the  spatial  directions  perpendicular  to  the  incident-wavefront 
propagation  direction,  the  OPD  profiles  are  in  this  case  functions  of  x  and  y  as  well  as  time, 
i.e. 

OPD(a;,  y,  t ) ,  (7) 

where,  for  simplicity  of  notation  in  what  follows  below,  x  and  y  are  the  spatial  coordinates 
normalized  by  the  large-scale  transverse  extent  L  of  the  flow  (cf.  figure  1),  and  OPD  will 
denote  the  optical  path  difference  normalized  by  (An)ref  L  with  (An)ref  denoting  the  reference 
(freestream-based)  refractive-index  difference.  To  be  able  to  vary  the  anisotropy  of  the  OPD, 
which  will  be  needed  in  §§  3  and  4,  we  propose  the  use  of  scaling  factors  (stretching  factors)  for 
both  the  OPD  and  the  spatial-extent  coordinates,  i.e.  so  that  the  anisotropic  transformation 
of  the  optical-wavefront  surface  can  be  represented  as 

OPD(x,y,t)  — *  aOPD  OPD(ax  x ,  ay  y ,  t) ,  (8) 

where  q;opd  is  a  scaling  factor  for  the  OPD,  and  {ax,ay}  are  scaling  factors  for  the  spatial- 
extent  coordinates  normal  to  the  incident-wavefront  propagation  direction.  These  scaling 
factors  can  each  be  greater  than  or  less  than  unity.  Two  examples  are  shown  in  figure  2 
where  different  values  of  q;opd  are  employed  for  the  same  wavefront  or  OPD  isosurface.  We 
note,  with  regard  to  equation  8,  that  only  two  wavefront-anisotropy  parameters  are  needed 
(only  the  ratios  matter)  and,  therefore,  we  define  a  (vector)  aerooptical-wavefront-anisotropy 
(AWA)  parameter  given  by 

_  f  <*opd  aopD  1  /nN 

in  terms  of  the  ratios  of  the  scaling  factor  of  the  OPD  to  each  of  of  the  scaling  factors  of 
the  transverse  spatial  coordinates.  For  a  wavefront  in  2-D  space,  e.g.  in  figures  3a  and  b, 
only  one  scaling  factor  is  needed,  e.g.  the  AWA  parameter  becomes  a  =  oopd/o^.  The 
proposed  AWA  parameter  a  enables  the  variation  of  the  effective  anisotropy  of  aerooptical 
wavefronts,  which  can  be  used  to  examine  their  structure  as  a  function  of  scale  as  will  be 
shown  in  §§  3  and  4.  We  note  that  the  AWA  parameter  can  be  used  not  only  in  the  study 
of  the  OPD  isosurfaces  but  also  in  the  study  of  the  OPD  field  itself,  if  needed. 

As  mentioned  above,  the  present  anisotropy  considerations  are  applicable  to  both  the 
wavefront  itself  and  to  the  OPL  or  OPD  profiles.  The  proposed  AWA  parameter,  therefore, 
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is  also  applicable  to  both.  In  practice,  it  is  important  to  appreciate  that  while  the  OPL 
isosurfaces  represent  the  wavefronts,  what  is  more  straightforward  to  measure  experimentally 
is  the  OPL  or  OPD  profile,  i.e.  the  variation  of  the  OPL  in  the  reference  measurement  plane 
of  a  detector  or  along  a  reference  line.  It  is  helpful,  however,  to  also  appreciate  that  the  OPL 
isosurfaces  are,  after  all,  derivable  from  the  OPL  field  assuming  the  OPL  field  is  known.  For 
example,  considering  wavefronts  that  emerge  from  a  turbulent  shear  flow  and  subsequently 
propagate  in  a  near-field  region  of  uniform  reference  refractive  index,  the  wavefronts  will 
propagate  at  a  constant  speed  in  the  uniform-index  region.  As  long  as  this  region  is  in 
the  near  field,  i.e.  without  having  to  consider  far-field  propagation  effects,  the  normalized 
OPL  profiles  are  geometrically  equivalent  to  the  OPL  isosurfaces.  The  wavefront- anisotropy 
considerations  above,  and  the  AWA  parameter,  are  equally  applicable  to  the  profiles  as  well 
as  isosurfaces  of  the  OPL  or  OPD  fields. 

3.  Anisotropic  box-counting  technique  for  characterization  of  aerooptical  wave- 
fronts 

To  investigate  the  physical  structure  of  optical  wavefronts,  a  useful  approach  is  to  con¬ 
duct  box  counting  which  enables  an  examination  of  the  manner  in  which  the  wavefronts’ 
geometrical  features  vary  with  scale  (e.g.  Catrakis  et  all7\  Catrakis  &  Aguirre18).  For  each 
wavefront,  a  bounding  box  or  rectangle  is  identified  and  partitioned  successively  into  smaller 
boxes.  The  number  of  boxes  containing  the  wavefront  is  counted  at  each  box  scale  A.  This 
coverage  count,  denoted  as  Nd( A)  in  d  dimensions,  can  be  used  to  quantify  the  geometric 
structure  as  a  function  of  scale.  At  the  largest  scale  Amax,  the  coverage  count  is  unity,  i.e. 
Nd{ Amax)  =  1.  We  sometimes  also  denote  the  largest  scale  as  J,  i.e.  the  bounding-box  scale. 
One  can  normalize  the  coverage  count  by  the  total  number  of  boxes  available  at  each  scale 
A  and  this  is  known  as  the  coverage  fraction, 

m)  =  (tM  Nd( a),  (io) 

where  0  <  Fd( A)  <  1  and  the  limiting  values  correspond  to  the  smallest  and  largest  scales, 
respectively.  The  physical  meaning  of  the  coverage  fraction  is  that  it  is  the  geometric  prob¬ 
ability  of  finding  a  part  of  the  wavefront  in  a  A-scale  box.  A  larger  value  of  Fd( A),  at  a  given 
scale,  indicates  a  more  irregular  wavefront  structure  at  that  scale. 
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Highly-anisotropic  behavior  can  be  expected  for  aerooptical  wavefronts  (cf.  §2)  and, 
because  of  this,  the  scale  A  must  be  treated  as  a  vector  in  general,  as  was  proposed  for 
turbulent  fluid  interfaces  in  general  by  Catrakis19.  In  the  present  case, 

A  =  {Aopd,  Ax,  Ay}  ,  (11) 

with  the  (scalar)  box  size  A  denoting  in  the  general  case  the  magnitude  of  A.  In  3-D  space, 
for  example,  we  can  write  this  as 


A  =  (Aopd  K  Ay)1/3  ,  (12) 

where  Aopd  denotes  the  extent  of  each  box  along  the  OPD  axis  and  {As,  Ay}  denote  the 
spatial  extents  of  each  box  in  the  {x,  y}  directions.  The  advantage  of  a  scale  magnitude 
defined  as  in  equation  12  is  that  A3  can  be  identified  as  a  volume  occupied  by  each  A-sized 
box. 


The  most  important  quantity  that  the  box-counting  technique  produces  is  the  coverage 
dimension  Dd( A)  as  a  function  of  the  scale  magnitude  A.  The  coverage  dimension  can  be 
computed  as  the  fractional  increase  in  the  coverage  count  per  unit  fractional  decrease  in 
scale,  or  equivalently  the  logarithmic  derivative  of  the  coverage  count  with  scale,  i.e. 


Dd(  A)  = 


dNdW/Nd{X) 

dA/A 


d  log  iVrf(A) 
dlog  A 


(13) 


and  can  be  interpreted  as  a  fractional  (fractal)  dimension  whose  departure  from  the  topolog¬ 
ical  dimension  dt  of  the  wavefronts  quantifies  how  complex  the  wavefront  structure  is.  The 
coverage  dimension  must  be  in  the  range 


dt  <  Dd{ A)  <  d,  (14) 

for  optical  wavefronts.  In  particular,  at  the  smallest  scales  and  largest  scales  respectively 
the  coverage  dimension  reaches  those  limits,  i.e. 

Drf(A)  *  d^  as  A  *  ^  Am,n ,  and ,  D^(A)  ^  d  as  A  ♦  Amax  *  (lb) 

For  optical  wavefronts  in  3-D,  as  indicated  in  figure  1,  these  limiting  values  become  <k  —  2 
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and  d  =  3  respectively.  For  wavefronts  in  2-D,  cf.  figure  3,  these  values  become  dt  =  l 
and  d  =  2  respectively.  It  is  also  helpful  to  think  of  the  relative  coverage  dimension  of  the 
wavefronts,  i.e.  the  coverage  dimension  relative  to  the  topological  dimension,  which  is  given 
by  Dd( A)  —  dt  in  general,  or  D2( A)  —  1  in  2-D  for  example.  This  is  useful  for  comparing 
wavefront  structure  in  different  dimensions,  e.g.  in  3-D  vs.  2-D,  because  the  values  of  the 
relative  coverage  dimension  must  range  from  0  to  1  independent  of  whether  the  wavefront  is 
in  a  3-D  or  2-D  region  of  the  flow,  cf.  equation  15.  Recalling  that  at  large  Reynolds  numbers 
one  may  expect  similarity  properties  over  a  wide  range  of  scales  (cf.  equation  3),  another 
important  use  of  the  coverage  dimension  is  that  it  can  reveal  such  similarity  properties, 
including  the  presence  of  self-similarity. 

If  the  coverage  dimension  Dd(X)  is  a  constant  over  a  wide  range  of  scales,  then  this  would 
mean  that  the  wavefront  physically  has  scale-independent  structure  in  that  range  of  scales 
and  then  the  wavefront  could  be  modeled  using  this  self-similarity.  Specifically,  for  self¬ 
similar  behavior  the  coverage  dimension  would  be  Dd(X)  =  Dd  =  const,  and  the  structure  of 
the  wavefront  could  be  extrapolated  to  larger  Reynolds  numbers  (cf.  equation  2)  using 

m a)  ~  ,  tie) 

over  a  range  of  scales  Amjn  <  A  <  Alt,  cf.  equation  3.  The  value  of  Dd  would  reflect  the 
complexity  of  the  wavefront  structure  that  statistically  repeats  itself  over  a  range  of  scales. 

If,  however,  the  coverage  dimension  is  found  to  be  a  continuous  function  of  scale  through¬ 
out  the  range  of  scales,  then  the  coverage  count  would  no  longer  be  a  power  law  but  instead 
would  behave  as 

U^max  dA'  1 

Dd(X')  ^  j  ,  (17) 

as  can  be  readily  seen  by  inverting  equation  10.  Such  scale-dependent  behavior,  where  for 
example  the  coverage  dimension  could  be  increasing  with  increasing  scale  over  a  range  of 
scales,  would  mean  that  the  complexity  of  the  wavefront  varies  with  scale.  This  would 
mean  that  more  complicated,  scale-dependent  models  would  be  needed  to  extrapolate  the 
small-scale  structure  of  optical  wavefronts  to  larger  Reynolds  numbers.  We  note,  however, 
that  scale-dependent  behavior  on  the  basis  of  Dd( A)  alone  may  not  necessarily  imply  lack 
of  self-similarity.  Recent  work  has  indicated  that  other  scale-local  measures  may  be  needed 
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to  discern  self-similar  behavior  in  those  cases  where  Dd{ A)  appears  to  be  scale  dependent 
(Catrakis,  Aguirre,  &  Ruiz-Plancarte5).  In  the  cases,  however,  where  there  is  clear  evidence 
of  scale  independence  in  Dd(X),  this  quantity  is  sufficient  to  detect  self-similarity. 

As  emphasized  in  §2,  aerooptical  wavefronts  can  be  expected  to  exhibit  a  strong  anisotropy. 
A  more  general  box-counting  approach  therefore  is  needed,  to  be  useful  for  aerooptics,  where 
one  can  examine  the  box-counting  behavior  of  the  wavefronts  allowing  for  the  possibility 
that  they  may  be  scale  independent,  i.e.  self-similar,  in  an  anisotropic  manner.  Self-similar 
anisotropic  behavior  is  sometimes  also  called  self-affine  behavior.  Irregular  self-affine  objects 
exhibit  structure  which  is  scale  independent  but  in  a  manner  that  depends  on  the  direction. 
In  other  words,  self-affine  objects  are  anisotropically  scale  independent  whereas  self-similar 
objects  are  isotropically  scale  independent. 

We  propose  an  anisotropic  extension  of  the  box-counting  technique  whereby  the  AWA 
parameter  a  of  §2  becomes  a  parameter  in  the  coverage  count  and  consequently  in  the 
coverage  fraction  as  well  as  in  the  coverage  dimension,  i.e. 

Nd( A;  a.) ,  and ,  Fd( A;  a) ,  with ,  Dd{ A;  a)  =  -  -  ~  •  (18) 

The  introduction  of  the  AWA  parameter  a  in  the  coverage  quantities  makes  possible  a 
practical  examination  of  anisotropic  structure.  We  note  that  an  equivalent  way  to  proceed 
would  be  to  treat  the  scale  A  as  a  vector  A  (cf.  equation  11)  and  then  to  define  the  coverage 
dimension  through  an  extension  of  equation  13.  This,  however,  would  lead  to  coverage- 
dimension  vectors  expressed  as  partial  derivatives  of  the  coverage  count.  The  proposed 
approach,  i.e.  using  the  quantities  in  equation  18,  offers  instead  a  much  simpler  way  to 
characterize  the  anisotropic  structure  of  aerooptical  wavefronts,  with  a  coverage  dimension 
that  is  scalar  valued  rather  than  vector  valued. 

As  an  example  in  2-D,  figures  4a  and  b  show  the  same  optical  wavefront  at  two  different 
values  of  the  AWA  parameter  a  with  the  boxes  necessary  to  cover  the  wavefront  in  each  case 
for  several  different  scales.  The  corresponding  coverage  results  are  shown  in  figures  5a,  b, 
and  c.  Note  that  the  choice  of  the  values  of  the  AWA  parameter,  in  figures  4a  and  b,  is  such 
that  they  both  have  the  same  scale  magnitudes,  i.e.  the  value  of  (oopl  ax)^2  is  the  same  for 
both  examples.  This  was  done  so  that  any  effects  of  the  largest  scales  are  removed  in  order  to 
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focus  on  the  small-scale  structures.  Significantly,  we  see  that  for  the  same  optical  wavefront 
examined  under  different  levels  of  anisotropy,  different  behavior  is  evident  at  the  small  scales 
in  figure  5c.  This  is  an  important  observation  because  it  illustrates  that  this  technique  is 
sensitive  to  the  presence  of  anisotropy  at  small  scales  and,  therefore,  one  can  expect  to  be 
able  to  use  it  to  examine  aerooptical  wavefronts  at  large  Reynolds  numbers  as  will  be  done 
in  §4.  It  is  also  important  to  appreciate  that  the  above  ideas  can  be  expected  to  be  useful 
to  examine  not  only  the  ensemble-averaged  behavior  of  the  instantaneous  wavefronts  but 
also  the  instantaneous  wavefront  structure  itself.  The  latter  is,  of  course,  crucial  in  practical 
aerooptics  applications. 

4.  Scale-local  structure  of  aerooptical  distortions  at  large  Reynolds  number  and 
high  compressibility 

In  this  section,  we  demonstrate  the  use  of  the  proposed  anisotropic  box-counting  tech¬ 
nique  for  the  characterization  of  the  scale-local  structure  of  aerooptical  distortions.  We  are 
particularly  interested  in  the  behavior  at  large  Reynolds  numbers  and  high  compressibility 
because  this  is  a  flow  regime  of  practical  interest  in  a  number  of  aerooptics  applications 
in  high-speed  flight,  and  also  because  previous  work  suggests  that  high-compressibility 
turbulent  flows  exhibit  important  differences  at  least  at  the  large  scales  when  compared 
to  incompressible  or  weakly-compressible  turbulent  flows  (e.g.  Samimy  &  Elliot20;  Pa- 
pamoschou21).  Less  is  known  about  the  effect  of  compressibility  on  the  small-scale  behavior  at 
large  Reynolds  numbers.  Computational  results  are  available  on  compressible  turbulent  flows 
(e.g.  Porter,  Woodward,  &;  Pouquet22)  but  are  presently  restricted  to  low  Reynolds  num¬ 
bers.  On  the  theoretical  side,  there  are  proposals  of  similarity  and  self-similarity  properties 
at  large  Reynolds  numbers  (e.g.  Sreenivasan4)  that  have  been  put  forth  for  incompressible 
turbulence.  However,  it  is  presently  not  clear  to  what  extent  such  incompressible-turbulence 
similarity  properties  are  applicable  at  compressible-flow  conditions  including  high  compress¬ 
ibility. 

To  examine  the  applicability  of  the  technique  proposed  in  §§  2-3,  we  employ  recent  wave- 
front  data  derived  from  experimental  studies  of  aerooptical  distortions  in  high-compressibility 
shear  layers  between  different  gases  (Dimotakis,  Catrakis,  &  Fourguette2).  A  discussion  of 
the  experimental  technique  can  be  found  in  the  companion  paper  (Catrakis  Sc  Aguirre14) 
with  a  more  complete  description  in  Dimotakis,  Catrakis,  &;  Fourguette2).  The  aerooptical 
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behavior  in  these  measurements  is  dominated  by  mixing  effects.  The  important  issue  of 
mixing  vs.  density  effects  in  aerooptics  is  discussed  in  Jumper  &  Fitzgerald3  and  also  in  the 
companion  paper  to  the  present  work  (Catrakis  &  Aguirre14).  For  the  data  employed  here, 
a  large-scale  Reynolds  number  of  Re  ~  106  and  a  convective  Mach  number  of  Mc  ~  1  were 
achieved  and  the  OPD  profiles  corresponding  to  wavefronts  emerging  normal  to  the  shear 
layers  were  derived  by  Fourier  filtering  the  refractive-index  field,  as  described  in  Dimotakis, 
Catrakis,  &  Fourguette2.  Examples  of  the  OPD  profiles  at  the  high-compressibility  flow  con¬ 
ditions  are  shown  in  figures  6a,  b,  and  c.  The  optical-wavefront  distortions  correspond  to 
the  aerooptical  interactions  in  a  two-dimensional  plane  transverse  to  the  shear  layer.  The 
incident  wavefronts  are  initially  linear,  i.e.  straight,  and  subsequently  propagate  transverse 
to  the  flow,  from  one  freestream  to  the  other,  through  two-dimensional  spatial  streamwise 
slices  of  the  flow.  The  wavefront  OPD  profiles  correspond  to  the  aerooptical  wavefronts  that 
emerge  from  the  turbulent  flow  region.  Although  the  wavefront  data  in  figures  6a  and  b  are 
not  fully  resolved,  the  resolution  is  sufficient  to  capture  a  large  fraction  of  the  full  range  of 
scales  at  these  large  Reynolds  numbers.  This  fraction  of  the  range  of  scales  is  estimated  to  be 
~  500  :  1  based  on  the  camera  pixel  resolution  which  was  employed  to  record  the  flow  images 
from  which  the  OPD  profiles  were  derived.  Keeping  in  mind  that  similarity  properties  are 
not  expected  throughout  the  entire  range  of  scales  smaller  than  the  large  scale,  but  rather 
in  a  fraction  of  those  scales  (cf.  equation  3),  these  data  can  be  expected  to  be  useful  to 
examine  the  presence  of  similarity  or  self-similarity  in  the  wavefront  structure. 

For  the  shear-layer  data,  the  OPD  profiles  are  functions  of  the  streamwise  coordinate  x 
normalized  by  the  transverse  extent  L  of  the  large-scale  shear  in  this  flow,  i.e. 

OPD  (|)  .  (19) 

The  incident  wavefronts  are  assumed  to  be  planar  (linear  in  the  planar  slice  of  the  flow 
studied)  so  that  the  coverage  dimension  is 

D2(  A)  =  1 ,  (20) 

prior  to  propagation  of  the  wavefronts  through  the  shear  layer.  Subsequent  to  propagation 
and  distortion  of  the  wavefronts,  the  relevant  AWA  parameter  associated  with  the  OPD 
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profile  is 


a  = 


^OPD 


(21) 


(cf.  §2  and  equations  8-9),  corresponding  to  the  two-parameter  stretched-OPD  profile,  i.e. 


c^opd  OPD 


(22) 


As  pointed  out  in  §2,  the  two  parameters  oopd  and  ax  are  equivalent  to  the  single  AWA 
parameter  a  =  (>opd /otx.  One  parameter  is  sufficient,  for  wavefronts  in  2-D  space,  to  permit 
an  examination  of  the  wavefront  structure  at  different  levels  of  anisotropy. 


Box  counting  was  performed  on  the  high-compressibility  shear-layer  wavefront  profiles  for 
several  different  values  of  the  AWA  parameter  a,  over  the  range  0.3  <  a  <  3,  i.e.  spanning 
one  order  of  magnitude  in  the  AWA  parameter.  For  each  value  of  a,  the  wavefront  was 
stretched  and  subsequently  partitioned  into  boxes  (tiles  in  2-D)  so  that  the  box  count  was 
computed  as  a  function  of  scale,  i.e. 

N2( A;  a),  (23) 

where  the  subscript  2  denotes  that  these  data  are  spatially  two-dimensional.  From  the 
coverage  count,  we  computed  the  coverage  fraction  and  the  coverage  dimension  as  a  function 
of  scale,  i.e. 

F2(A;a)  =  (-M  N2( A; a),  and,  D2{ A; a)  =  _  dlo^lV^a)  (24) 
\  ^max  /  d  log  A 

Two  important  observations  were  made: 


(a)  the  small-scale  wavefront  behavior  depends  strongly  on  the  AWA  parameter  a,  and 

(b)  there  exists  a  particular  value  of  the  AWA  parameter  a  =  a*  for  which  the  small-scale 

wavefront  structure  is  found  to  be  scale  independent  over  a  range  of  scales. 

The  first  observation  confirms  the  idea  in  §  2  that  aerooptical  wavefronts  are  highly  anisotropic 
and  shows  that  the  technique  described  in  §  3  is  capable  of  detecting  the  dependence  of  the 
small-scale  structure  on  this  anisotropy.  The  second  observation,  i.e.  that  there  is  scale 
independence  over  a  range  of  scales  for  a  particular  value  a*,  which  is  found  to  be  a*  ~  0.9 
for  the  present  data,  is  particularly  important.  The  results  are  shown  in  figures  7a,  b,  and 
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c,  that  depict  respectively  the  coverage  count  N2( A;  a*),  coverage  fraction  F2( A;  a*),  and  the 
relative  coverage  dimension  D2(X;a*)  -  1  for  the  critical  value  of  the  AWA  parameter.  For 
values  of  a  smaller  than  or  greater  than  a*,  qualitatively-different  behavior  was  observed. 
Specifically, 

(a)  for  a  <  a*,  the  wavefront  structure  is  scale  dependent  with  relatively-smaller  values  of 

the  coverage  dimension  at  intermediate  scales, 

r 

(b)  there  exists  a  particular  value  of  the  AWA  parameter  a  =  a*  for  which  the  small-scale 

wavefront  structure  is  found  to  be  scale  independent  over  a  range  of  intermediate  scales. 

(c)  for  a  >  a*,  the  wavefront  structure  is  scale  dependent  with  relatively-larger  values  of 

the  coverage  dimension  at  intermediate  scales. 

The  significance  of  the  observed  behavior  at  a  =  a*  is  that  this  shows  the  presence  of  an 
intermediate  range  of  small  scales  in  which  the  wavefront  is  self-similar,  i.e.  a  range  of  scales 
in  which  the  coverage  dimension  is  nearly  a  constant: 

D2{\\a*)  ~  const.  ~  1.27.  (25) 

The  uncertainty  in  the  coverage  dimension  is  estimated  to  be  ±0.025,  estimated  from  the 
variations  between  individual  realizations.  The  D2  ~  1.27  value  is  found  for  the  present  flow 
conditions  of  high  compressibility  and  large  Reynolds  number.  At  other  flow  conditions,  e.g. 
different  compressibility  or  Reynolds  numbers,  it  may  exhibit  a  dependence  but  this  would 
be  expected  to  be  a  weak  dependence  as  long  as  the  flow  is  fully  developed.  It  can  be  noted 
that  the  D2  ~  1.27  value  is  lower  than  the  D2  ~  1.35  value  reported  for  incompressible 
flows  (e.g.  Sreenivasan4),  and  this  is  consistent  with  expectations  that  the  dimension  should 
decrease  with  increasing  compressibility  (e.g.  Smits  &  Dussauge12).  The  present  finding  can 
be  used  to  extrapolate  the  behavior  to  even  larger  Reynolds  numbers  at  high  compressibility, 
at  least  in  shear  layers,  using 

~  (^y  27,  (26) 

(cf.  equation  16).  From  figure  7c,  we  see  that  the  self-similar  behavior  is  evident  in  the 
measurements  over  a  range  of  scales  that  span  approximately  one  decade.  Even  though  the 
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Reynolds  number  associated  with  the  shear  layer  is  large,  i.e.  Re  ~  106,  and  corresponds  to 
a  full  range  of  flow  scales  given  by  Amax/Amin  =  Re 3//4  ~  3  x  104,  the  range  of  scales  over 
which  similarity  is  expected  is  smaller.  For  example,  a  conservative  estimate  (Dimotakis6) 
is  the  ratio  of  the  Liepmann-Taylor  scale  Alt  to  the  Kolmogorov  scale  Ak  which  grows 
slowly  with  Reynolds  number,  specifically  Alt/Ar  r\_/  Re1/4  n-i  10,  cf.  equations  2  and  3 
in  §2.  This  estimate,  at  least  of  the  extent  of  the  range  of  self-similar  scales,  is  closer  to 
the  presently-observed  behavior  (cf.  figure  7c).  Higher-resolution  measurements  would  be 
needed  to  discern  this  issue  further,  i.e.  whether  the  range  of  scales  exhibiting  self-similar 
behavior  extends  to  smaller  scales  yet.  We  note  that  the  critical  value  a*  of  the  AWA 
parameter  might  depend  to  some  extent  on  the  Reynolds  number  and/or  Mach  number, 
and  possibly  on  the  flow  boundary  conditions.  However,  as  long  as  the  Reynolds  number  is 
large,  i.e.  meets  the  criterion  for  fully-developed  turbulence  (cf.  §2),  one  can  expect  that 
the  small-scale  structure  of  the  wavefronts  and  the  critical  value  of  the  AWA  parameter  will 
be  at  most  weakly  dependent  on  these  flow  parameters.  This  is  because  the  small-scale 
structure  at  large  Reynolds  numbers  is  at  most  weakly  sensitive  to  the  large-scale  behavior 
(e.g.  Sreenivasan4,  Dimotakis6,  and  references  therein). 

The  present  finding  of  scale  independence  of  the  small-scale  optical-wavefront  struc¬ 
ture  over  a  range  of  scales  at  high  compressibility  provides  experimental  evidence  that 
self-similarity  is  present  in  aerooptical  interactions  at  high-compressibility  flow  conditions. 
Computational  studies  of  compressible  flows  at  low  Reynolds  numbers  have  indicated  inertial- 
range  scaling  based  on  power  spectra  (e.g.  Porter,  Woodward,  &  Pouquet22)  and  this 
supports  the  present  finding  of  self-similarity,  although  spectral  behavior  alone  cannot  be  re¬ 
lated  uniquely  to  the  physical  structure  of  the  fluid  interfaces.  The  present  observation 
of  self-similar  behavior  in  the  physical  structure  of  the  distorted  aerooptical  wavefronts 
indicates  that  turbulent  fluid  interfaces  at  high  compressibility  also  possess  self-similarity 
properties.  Self-similar  behavior  of  the  fluid  interfaces  would  be  sufficient  to  explain  the  ob¬ 
served  self-similar  aerooptical  behavior,  with  anisotropy  inherent  in  the  wavefront  OPD,  cf. 
§2.  Quantification  and  modeling  of  the  physical  structure  of  optical  wavefronts,  in  real  space 
as  opposed  to  Fourier  space,  requires  tools  such  as  the  technique  presently  demonstrated. 
The  present  results  demonstrate  that  optical-wavefront  self-similarity  can  be  detected  with 
the  proposed  approach  using  the  AWA  parameter. 
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5.  Conclusions  and  implications 

The  proposed  method  of  anisotropic  box  counting  provides  a  technique  useful  for  the 
scale-local  characterization  of  the  physical  structure  of  aerooptical  wavefronts,  i.e.  as  a 
function  of  scale.  The  aerooptical-wavefront-anisotropy  (AWA)  parameter,  defined  as  a 
ratio  of  scaling  factors  for  the  wavefront  distortions  and  for  the  transverse  extent  of  the 
wavefronts,  combined  with  an  anisotropic  generalization  of  box  counting  enables  a  scale- 
local  examination  of  wavefronts  at  varying  degrees  of  anisotropy.  This  is  particularly  useful 
for  the  study  of  aerooptical  wavefronts  at  large  Reynolds  numbers.  Demonstration  of  this 
technique  to  experimental  data  in  large-Reynolds-number  (Re  ~  106)  high-compressibility 
(Mc  ~  1)  turbulent  shear  layers  shows  the  presence  of  anisotropic  self-similarity,  or  self- 
affinity,  in  the  wavefront  structure  over  a  range  of  small  scales  spanning  nearly  a  decade, 
as  detected  by  varying  the  AWA  parameter  and  conducting  a  scale-local  examination  of  the 
distortions.  Even  though  the  data  employed  in  this  study  are  dominated  by  mixing  effects 
as  opposed  to  solely  density  effects,  the  anisotropy  considerations  and  the  AWA  parameter 
can  be  expected  to  be  useful  in  both  cases. 

The  present  finding  and  technique  have  helpful  consequences  for  large-Reynolds-number 
experimental  studies  as  well  as  modeling  and  computational  studies.  The  data  show  the 
utility  of  the  technique  to  detect  scaling  in  large-Reynolds-number  flow  experiments  using 
the  optical  behavior,  and  as  such  suggest  an  optical  means  to  probe  aspects  of  the  mixing 
behavior  at  large  Reynolds  numbers.  The  present  finding  also  provides  a  means  to  extrapo¬ 
late  the  small-scale  structure  of  optical  distortions  at  high  compressibility  to  larger  Reynolds 
numbers.  This  can  be  expected  to  be  useful  in  efforts  to  develop  computational  models  and 
simulations  of  optical  behavior  in  flows  at  large  Reynolds  numbers  and  high  compressibil¬ 
ity  (e.g.  Jones  &  Bender23).  The  present  self-similarity  model  of  the  wavefront  distortions 
suggests  a  key  ingredient  to  help  construct  effective  sub-grid-scale  models  in  large-eddy  sim¬ 
ulations  of  aerooptical  interactions.  Small-scale  properties  of  high-compressibility  flows  are 
also  valuable  because  they  can  be  coupled  with  descriptions  that  address  the  large-scale 
behavior  of  high-compressibility  flows  (e.g.  the  companion  paper  by  Catrakis  &  Aguirre14), 
in  order  to  develop  frameworks  that  are  capable  of  describing  the  behavior  across  the  en¬ 
tire  range  of  scales.  Furthermore,  this  can  provide  the  means  to  compare  the  cumulative 
effects  of  distortions  at  small  scales  to  the  effects  at  the  large  scales.  These  implications  are 
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especially  important  in  applications  that  require  high-resolution  optical  imaging  and  beam 
propagation  in  high-speed  flight. 
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Figure  1.  Schematic  of  a  planar  optical  wavefront  (top)  incident  to  a  three-dimensional 
refractive-index  field  in  a  large- Reynolds-number  turbulent  compressible  flow  (middle),  and 
the  emerging  highly-irregular  aerooptical  wavefront  (bottom). 

Figures  2a  and  b.  Illustration  of  varying  the  effective  anisotropy  of  an  optical  wavefront  in  3- 
D  using  the  proposed  aerooptical-wavefront-anisotropy  (AWA)  parameter  a,  cf.  equations  8 
and  9.  The  same  wavefront  is  shown  with  |o:|  3>  1  (top)  and  with  |a|  C  1  (bottom). 

Figures  3a  and  b.  Example  of  an  optical  wavefront  in  2-D.  The  same  wavefront  is  shown  at 
two  different  values  of  the  AWA  parameter,  which  is  a  scalar-valued  parameter  in  2-D  and 
is  a  =  qjopd /cKx-  In  this  example,  the  large-scale  size  of  the  wavefront  bounding  box  is  the 
same  for  each  case,  i.e.  (ciopd  o^)1/2  is  the  same. 

Figure  4a.  Illustration  of  the  anisotropic  box-counting  technique  for  the  wavefront  as  shown 
in  figure  3a.  Coverage  boxes  are  shown  at  different  scales.  The  coverage  boxes  are  the  boxes 
needed  to  completely  contain  the  wavefront  at  each  scale. 

Figure  4b.  Illustration  of  the  anisotropic  box-counting  technique  for  the  wavefront  as  shown 
in  figure  3b. 

Figures  5a,  b,  and  c.  Comparison  of  the  coverage  count,  coverage  fraction,  and  coverage 
dimension  for  the  wavefront  as  shown  in  figures  3a  and  b.  The  dashed  curves  correspond  to 
figure  3a,  and  the  solid  curves  correspond  to  figure  3b.  The  differences  at  the  small  scales 
indicate  that  the  small-scale  behavior  is  sensitive  to  the  choice  of  the  AWA  parameter. 

Figures  6a,  b,  and  c.  Examples  of  experimental  wavefront  OPD  profiles,  OPD(x),  derived 
from  two-dimensional  refractive-index  field  measurements  in  large-Reynolds-number  (Re  ~ 
106)  high-compressibility  (Mc  ~  1)  shear  layers.2,14  The  OPD  and  spatial  coordinate  x 
denote  normalized  values.  The  normalization  employs  the  transverse  spatial  extent  of  the 
shear  layer  and  the  ambient  (freestream)  refractive-index  difference.  The  incident  wavefronts 
are  initially  linear,  i.e.  straight,  and  propagate  transverse  to  the  flow,  from  one  freestream 
to  the  other,  through  two-dimensional  spatial  streamwise  slices  of  the  flow.  The  wavefront 
OPD  profiles  correspond  to  the  aerooptical  wavefronts  that  emerge  from  the  turbulent  flow 
region. 
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Figures  7a,  b,  and  c.  Ensemble-averaged  results  for  the  instantaneous  structure  of  wavefront 
OPD  profiles  in  large-Reynolds-number  (Re  ~  106)  high-compressibility  (Mc  ~  1)  shear  lay¬ 
ers,  cf.  figures  6a,  b,  and  c:  (a)  coverage  function  A^(A;  a*),  (b)  coverage  fraction  A;  a*), 
and  (c)  relative  coverage  dimension  Z?2(A;  a*)  —  1.  The  value  of  the  critical  AWA  parameter 
a*  for  the  results  shown  is  a*  ~  0.9.  The  plateau  in  (c)  indicates  a  clear  presence  of  self¬ 
similarity  of  the  anisotropic  wavefront  structure  over  a  range  of  small  scales  spanning  nearly 
a  decade. 
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AIAA  Journal,  Catrakis,  Aguirre,  &  Nathman,  Figures  3a  and  b 
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